Localization in disordered superconducting wires with broken spin-rotation symmetry 
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Localization and derealization of non-interacting quasiparticle states in a superconducting wire 
are reconsidered, for the cases in which spin-rotation symmetry is absent, and time-reversal symme- 
try is either broken or unbroken; these are referred to as symmetry classes BD and Dili, respectively. 
We show that, if a continuum limit is taken to obtain a Fokker-Planck (FP) equation for the transfer 
matrix, as in some previous work, then when there are more than two scattering channels, all terms 
that break a certain symmetry are lost. It was already known that the resulting FP equation ex- 
hibits critical behavior. The additional symmetry is not required by the definition of the symmetry 
classes; terms that break it arise from non-Gaussian probability distributions, and may be kept in 
a generalized FP equation. We show that they lead to localization in a long wire. When the wire 
has more than two scattering channels, these terms are irrelevant at the short distance (diffusive 
or ballistic) fixed point, but as they are relevant at the long-distance critical fixed point, they are 
termed dangerously irrelevant. We confirm the results in a supersymmetry approach for class BD, 
where the additional terms correspond to jumps between the two components of the sigma model 
target space. We consider the effect of random n fluxes, which prevent the system localizing. We 
show that in one dimension the transitions in these two symmetry classes, and also those in the 
three chiral symmetry classes, all lie in the same universality class. 



I. INTRODUCTION 



Localization of single particle states in disordered 
fermionic systems has been a key concept underlying our 
understanding of disordered normal metals. The nature 
of the states depends crucially on the presence or absence 
of time reversal (TR) and spin rotation (SR) symmetries. 
Based on these considerations, three distinct Wigner- 
Dyson symmetry classes were identified a long time ago 
— orthogonal, unitary and symplectic (the names come 
from the random matrix theory of Wigner and Dyson 1 ). 

Recently, it has been realized that localization ideas 
are useful in describing a wider class of systems including 
superconductors and superfluids, various quantum Hall 
states and random-bond Ising models. A crucial part of 
the emerging picture is again the role played by symme- 
tries of the disordered single-particle Hamiltonians de- 
scribing these systems. In addition to the TR and SR 
symmetries, these new systems may possess additional 
symmetries that have to be taken into account. One such 
category is comprised of systems defined on bipartite lat- 
tices with the sublattice (also called chiral) symmetry. 
Another category, on which this paper is focussed, in- 
cludes Bogoliubov-de Gennes (BdG) Hamiltonians which 
describe quasiparticles in disordered superconductors at 
a mean field level. Although particle number is not 
conserved in the mean field description, the problem of 
finding the quasiparticle spectrum can be transformed 
to a particle-number conserving form, so that it reduces 
to solving a single-particle Hamiltonian that has certain 
symmetries. The spectra of Hamiltonians in all of these 
classes have an exact reflection symmetry around zero 
energy, and eigenstates near zero energy may have local- 
ization properties different from those of eigenstates in 



other parts of the spectra. 

A symmetry classification of disordered single-particle 
Hamiltonians has been established by Altland and Zirn- 
bauer (AZ) 2 ^. The symmetry class and the dimension- 
ality of a particular disordered system in principle de- 
termine which phases (localized or extended) and phase 
transitions may occur in its phase diagram. After ensem- 
bles of random matrices, which are the zero-dimensional 
cases, the next simplest situation arises in one dimension 
(ID), where various analytical methods are available. In 
the standard Wigner-Dyson classes all the single-particle 
states are localized in ID and this leads to an exponential 
decrease of the conductance of a ID wire with its length, 
that is, the so-called strong localization. 

Brouwer et al. in Ref. studied localization in ID (su- 
perconducting wires) using the Fokker-Planck (FP) scat- 
tering approach^. They considered all four BdG sym- 
metry classes (denoted C, CI, D, and Dili in Ref. |3), and 
obtained the corresponding four FP equations describing 
transport at zero energy. While for the two classes C 
and CI that have unbroken SR invariance the quasipar- 
ticle (thermal) conductance Q decays exponentially with 
the length L of the wire for large L, the situation was 
found to be quite different in classes D and Dili in which 
the SR invariance is absent. In the latter, the mean (Q) 
decays only algebraically to zero for large L and h\Q is 
not self-averaging, indicating a very broad distribution of 
the conductance and the absence of the exponential lo- 
calization of the quasiparticle states. This behavior was 
termed critical in Ref. |j. 

This surprising result was questioned by Motrunich et 
al. in Ref. [jj In this paper some models in symmetry 
classes D and Dili were studied using real-space renor- 
malization group methods, and it was argued that, for 
a generic distribution of disorder, at zero energy there 
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are localized phases separated by a phase boundary at 
which the localization length diverges. It was found that 
close to the critical point, the low-energy density of states 
(DOS) per unit length, v{e) 1 has a divergent power-law 
behavior with a non-universal varying exponent: 

Ke)oce-*, (1) 

as e — > + , with S < 1. At the critical point, one has 
<5 = 1, and there are logarithmic corrections of the Dyson 
forrrA 

v{e) oc earner 3 . (2) 

In addition, as parameters vary, it is possible for the 
phase boundaries to collide, in which case there is a mul- 
ticritical point with different critical properties. 

In Ref. 9, the FP scattering approach was used to study 
the DOS in disordered wires in all symmetry classes. Per- 
haps unsurprisingly, for the wires in classes D and Dili 
the result had a singularity of the critical form in Eq. J2J. 
In Ref. 10, Brouwer et al. tried to account for the discrep- 
ancies between the results of the two groups in terms of a 
crossover from the unitary class to class D within a par- 
ticular model. We believe that this account leaves the 
role of symmetries and universality in the problem un- 
clear, and that our work described here greatly clarifies 
these questions. 

In this paper we reconsider localization in supercon- 
ducting wires with broken spin-rotation symmetry, in 
classes BD and Dili in the AZ classification [we prefer 
the term BD in place of class D for more than zero dimen- 
sions, as the distinction between odd (B) and even (D) 
size Hamiltonian matrices disappears for infinite matri- 
ces], for a finite number N of propagating modes in the 
wires, and resolve the major issues. The central point 
concerns symmetries of the problem. The most popular 
type of models for disordered wires are those in which 
propagation along the wire is continuous, and this ap- 
plies in particular for weak disorder. We show that for 
the cases of classes BD and Dili, there are certain sym- 
metry operations on the spaces of transfer matrices (or on 
the probability distributions for the disorder) which have 
been overlooked previously. The conditions that define 
classes BD and Dili do not impose that the distribution 
be invariant under these operations. However, when one 
takes the standard limit that results in the FP equation, 
the leading symmetry-breaking term is a differential op- 
erator of order N, and terms of order greater than two 
disappear in the limit (these are related to deviations 
of the disorder distribution from Gaussian). That is, the 
FP equations for classes BD and Dili studied by Brouwer 
et a&& possess more symmetry than the general class of 
distributions does. For N = 1, the leading symmetry- 
breaking term is a drift (first-order) term which produces 
localization if it is present. We show that for a long wire, 
the N channel case reduces to the N = 1 situation, and 
the leading symmetry-breaking term reduces to the drift 
term. Then generically, breaking the additional symme- 
try produces localization. In a generalization of the FP 



equation that keeps higher-order operators, operators of 
order higher than 2 would be irrelevant at short distance 
[in renormalization group (RG) language]. Since they 
are nonetheless relevant at large distances (they take the 
system off criticality) , they cannot be dropped, and so 
they are analogous to "dangerously irrelevant" operators 
in critical phenomena. In practice, to represent a system 
that will have both a finite number N of channels, and 
a finite mean free path as an ultraviolet cutoff, the FP 
limit is not taken, and symmetry-breaking terms must 
be retained. The behavior near the critical point is uni- 
versal, and is governed by the N = 1 case. The results 
agree with those of Motrunich et ai£. 

In addition to demonstrating these points with the FP 
approach, we also give the supersymmctric representa- 
tion of the wire for class BD, and show how the phenom- 
ena appear there. The wire is represented by a finite-A 
analog of the supersymmetric nonlinear sigma model (in 
ID), which is known to have a target space that consists 
of two connected components in the cases of classes BD 
and Dili 11 . We show that the terms in the distribution 
of disorder that break the additional symmetry produce 
a jump from one component to the other whenever they 
occur in the expansion of a disorder average. These are 
domain walls in the sigma model; topologically, in ID 
these domain walls are simply points. The supersymme- 
try approach is complementary to the FP equation; it 
produces results for the moments of Green's functions, 
and is used here to study the DOS. 

For the more general situation in which the transfer 
matrix does not have to evolve continuously, as in a lat- 
tice (or network) model, it is possible for an amplitude to 
pick up a factor of —1 during a time step, which would 
correspond to encircling a flux of it in a lattice model; 
such an effect is consistent with the conditions that define 
classes BD and DHL It was already shown— in the case of 
class BD with N = 1 that if such events are distributed 
independently and uniformly along the wire, the drift 
term in the FP equation is effectively suppressed, and 
the system is again critical. We show here that this also 
occurs for the case of general A, and then on scales larger 
than the separation of these events, the A-channel prob- 
lem effectively maps onto the analysis as given in Ref.S 
and is always critical even when the symmetry-breaking 
terms are present. In the supersymmetry approach, do- 
main walls are suppressed in this situation. Thus domain 
walls are essential to produce localization for all A, as is 
found also in two dimensions for class BDiA*i£. 

It is striking that the universal results for the sym- 
metry classes BD and Dili turn out to be the same, and 
also are the same as those for the chiral classes. We show 
that this "super-universality" can be understood, as the 
A = 1 channel cases all possess hopping Hamiltonian 
models that can be mapped onto each other exactly. 

The rest of the paper is structured as follows. In Sec- 
tion ^ we provide a brief account of the scattering ap- 
proach to coherent transport in wires, stressing its sym- 
metry aspects, and including the role of non-Gaussian 
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distributions, especially in the case of classes BD and 
DHL In Section 11111 we review the standard FP equa- 
tion resulting from this approach, and its solution. We 
generalize the FP equation to include the leading terms 
that break the extra symmetries. We simplify the result- 
ing generalized FP equation in the asymptotic regime of 
long wires and solve it, obtaining the scaling function of 
the logarithm of the conductance near the critical point, 
and exhibiting two localization lengths, the typical local- 
ization length for the typical values (and moments) of 
the logarithm of the conductance, and the mean local- 
ization length for the moments of the conductance itself. 
In Section llVl we describe aspects of the supersymmetry 
method for wires in the BD class. This gives us the scal- 
ing of the mean localization length and the mean density 
of states near the critical point. In Section we explain 
what occurs in more general models in which the propa- 
gation along the wire is discontinuous, allowing fluxes of 
it to be randomly inserted. The mappings that establish 
super-universality are given in Section IVII Some final 
discussion of universality is given in the conclusion, Sec- 
tionlVni Finally, many technical details are presented in 
Appendices. 



II. SYMMETRY ASPECTS OF THE TRANSFER 
MATRIX 

In this section we consider the symmetry aspects 
within the transfer matrix formulation, which will be 
used in both the FP and supersymmetry approaches, 
and show that for classes BD and Dili with more than 
two scattering channels there are additional symmetries 
present in a Gaussian distribution of disorder that are 
not generically there. In order to make the key points in 
our analysis, it will be useful to reconsider briefly the ba- 
sics of the transfer matrix approach to one-dimensional 
(ID) quantum transport (for more details, see Ref. [13J). 

It is convenient to pass from the single-particle Hamil- 
tonian TL to a scattering description, by considering 
waves of a given energy that enter the system at either 
end, are scattered and then emerge at the same or the 
opposite end. We will assume that only a finite number 
of "channels" are relevant at the energy considered. For 
the non-standard ensembles, the energy chosen has the 
special value e = (this corresponds to calculating the 
transport properties of the ground state of the supercon- 
ductor). The scattering can then be represented by an S- 
matrix that maps incoming to outgoing waves, or, better, 
by a transfer matrix M' that maps the in- and outgoing 
waves at one end (say, the right) to those at the other 
(the left); both of these matrices are 2A' x 2A'. The ad- 
vantage of the transfer matrix is that ID systems can be 
composed end-to-end by joining together short elements, 
and multiplying the respective transfer matrices, see Fig. 
^ Taking the latter as random, identically-distributed, 
statistically-independent matrices, leads to an equation 
that describes the evolution of the probability distribu- 




FIG. 1: The transfer matrix for the wire is composed of m 
steps, with independently-distributed increments Tj, j = 1, 
. . . , m. 

tion for the transfer matrix with increasing length of the 
system. 

The transfer matrix is related to the Landauer conduc- 
tance of the system. The Landauer-Buttiker formula (see 
Ref. [lj for history) for the two-probe (electrical) conduc- 
tance reads 

Q = ixttf, (3) 

where t is the A' x AT' transmission matrix for waves 
entering at one end to exit at the other. In supercon- 
ductors, charge (or quasiparticle number) is not a good 
quantum number at mean field level, but the Landauer 
conductance is still relevant for the thermal conductance, 
since heat is transported by quasiparticles: the thermal 
conductance is given by ^-k B TQ. The conductance can 
be rewritten in terms of the transfer matrix as 

Q = 2tr (2 + M'M' ] + (M'M' ] )- X )~ X . (4) 

Thus it is a trace of a function only of M'M (or of 
Af'tM'). 

The mathematical properties of M' are as follows. 
Because of conservation of probability flux ("current") 
in time-independent scattering, the S'-matrix is unitary, 
which implies that the transfer matrix is pseudounitary: 

M 1 ^ KM' = A, (5) 

where A = ^ ^ j- ^ . [The first A' indices refer 

to waves entering at the right (or leaving at the left), 
the second to those leaving at the right (or entering at 
the left).] Thus M' is a member of the noncompact Lie 
group U(A',A'). For the unitary ensemble, this is the 
only condition, but for the other ensembles 2 , additional 
conditions are imposed on M' . They result from the 
symmetries of the Hamiltonian TL that define the sym- 
metry class, and have the effect of restricting M' to a 
noncompact subgroup G of U(A', A'). In some cases, 
the representation of G in U(A', A') may be reducible 4 . 
That is, there exists a unitary transformation to a basis 
in which the full 2A' x 2N' transfer matrix M' becomes 
block diagonal, each of the d blocks being (for symme- 
try reasons) the same 2A x 2A transfer matrix M (or in 
some cases, its adjoint, or transpose), so A' = dN, where 
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d depends on the symmetry class. We will call N (rather 
than dN) the number of channels, as it is the number of 
independent variables propagating in each direction. 

In this section, and in the remainder of this paper ex- 
cept Section we aim to model the situation in which 
the propagation along the wire is continuous. For cal- 
culational convenience, we nonetheless represent this by 
discrete steps, as multiplication of transfer matrices for 
each step. The underlying continuity implies that after 
each step, M' lies in G and can be continuously con- 
nected to the identity, and thus G must be a connected 
group. For class BD, M' — M can be taken as real, 
so G is SO (A, A) with N' = N (as always, the prefix 
S means unit determinant, while the subscript denotes 
the component containing the identity). For class Dili, G 
is isomorphic to SO (2 A, C) (N' = 2N). For the purposes 
of most of this paper, these statements can be taken as 
the definitions of these symmetry classes. Tight-binding 
models of superconductors, though maybe not strictly 
continuous, lead to these systems with N even in both 
symmetry classes: Pairing of spinless fermions (without 
time-reversal symmetry) gives class BD, while pairing of 
spin- 1/2 fermions with time-reversal symmetry, but with 
broken spin-rotation symmetry, gives class DHL Pairing 
of fermions with spin 1/2, but with spin-rotation and 
time-reversal symmetries broken, gives class BD with N 
divisible by 4. Our conventions differ from those of Ref.|j, 
where their N is half ours for class Dili, and one-quarter 
of ours for class BD because they considered the case of 
spin- 1/2 fermions. Similarly, some of the models studied 
in Ref. are in class BD with N = 2 or 4, and others 
are in class Dili with N = 2. Our analysis applies to all 
positive integer values of A in both classes, and one odd 
value, N = 1, will play a special role in the analysis. 

Mathematically, class Dili Hamiltonians for a ID wire 
can also lead to transfer matrices in SO(2A + 1, C). Here 
we would have N' = 27V+1 instead of N' = 2N. The be- 
havior of these systems is different from that of those we 
analyze here, and somewhat similar to the odd number of 
channels case of the symplectic (All) symmetry classi^; 
namely the conductance of a long wire approaches Q = 1, 
an effect that was found in Ref. [H These "N +1/2- 
channel" class Dili systems can be realized as edge states 
of certain class Dili superconductors in two space di- 
mensions. We do not consider these cases further in this 
paper. 

The formula Q for the conductance becomes, in terms 
of M, 

g = 2dti (2 + MM f + (MM*)" 1 )" 1 , (6) 

and is invariant under multiplication of M on the left 
or right by unitary matrices. Those that lie in G rep- 
resent some redundant aspect of the evolution of M, in 
the sense that some scattering does not affect the con- 
ductance. We can find a maximal unitary (hence com- 
pact) subgroup K of G, such that Q is invariant under 
M — > kiMkz, k\,k2 £ K. This suggests parameter- 
izing M so as to factor off this subgroup. For exam- 



ple, for class BD, K = SO (A) x SO (TV), and for Dili, 
K = SO (2 AT). Then any transfer matrix M for class BD 
can be written schematically as M — k\ak2, or explicitly 

M _(Vi 0\ /cosh A sinhA\ [V 3 0\ m 
V° V 2 J \smhX coshxj \0 V 4 ) ' 1 ' 

where V\, . . . , V4 € SO (AT). For class Dili, we have in- 
stead 

and Vi, V2 are in SO (2 A). Thus the group can be de- 
composed as G = KAK (the Cartan decomposition of 
G), where A stands for the matrices containing X in the 
above expressions. In both cases, X = diag(xi, . . . ,xn) 
is real and diagonal. These parametrizations are also 
called radial or polar decompositions. The real param- 
eters Xi are called radial coordinates, and the matrices 
V\, . . . , V4, (resp., V\, V2) represent angular coordinates. 
Similar decompositions exist for other groups. A famil- 
iar example of the Cartan decomposition in a compact 
group is to take G = SU(2), K = U(l) (generated by the 
a z Pauli matrix), and then the three parameters in the 
decomposition are exactly the Euler angles that describe 
a rotation. 

If we identify elements of G modulo right multiplica- 
tion by elements of K, we obtain the coset space G/K. If 
we also identify points of this space modulo the action of 
K given by left multiplication, then we obtain the dou- 
ble coset space K\G/K. The X4 parametrize this double 
coset space K\G/K, and the Landauer conductance Q is 
a function only of these radial coordinates: 

N 1 

^ cosh 2 Xi ^ ^ 

To ensure that each point of K\G/K is labeled by a sin- 
gle set of Xi , there should be restrictions on the possible 
values of x;. We notice that it is possible to map the 
set A into itself by conjugation by certain elements of K, 
a — > kak' 1 . Since left and right multiplication by K is 
to be divided out, these give the same point in the dou- 
ble coset space. The group of such operations is called 
the Weyl group of the Cartan decomposition in question. 
By inspecting K for the cases of class BD and Dili, we 
see that it is possible in this way to permute the XiS, 
and to reverse the sign of an even number of XiS [because 
for class BD, K involves SO (AT), not O(N), and simi- 
larly for Dili]. Consequently, a unique parametrization 
for K\G/K is obtained by using the following domain 
(known as the Weyl chamber) for x^. 

C = {|sgi| <x 2 < ... <x N }, (10) 

where X\ can be either positive or negative. 

We want to note here that all the previous treatments 
of classes BD and Dili in ID based on the scattering 
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approach, including Refs. IH9 HT0I incorrectly assumed a 
different range for Xf 

< xi < x 2 < ■ ■ ■ < x N . (11) 

This range is the one appropriate for all the standard 
(Wigner-Dyson) symmetry classes, and for the BdG 
classes C and CI in the AZ classification. For the chi- 
ral classes, the range is 

X\ < X2 < ■ ■ ■ < xn- (12) 

The property that at least one of the radial coordinates 
(xi) can take negative values is shared by the classes BD, 
Dili, and the chiral ones. Ultimately, as we show below, 
this allows phase transitions to occur in these classes. 
Within the FP approach, a transition at which the lo- 
calization length diverges occurs when the mean of an Xi 
(in the limit of a long wire) passes through zero. There 
is one such transition in classes BD and Dili, but N 
in the chiral classes. (Related distinctions between the 
two phases, imply ing that there must be a transition, ap- 
peared in Refs. I7lll7l) The distinction in the ranges for 
Xi is crucial to our symmetry analysis also, as we will see 
that the additional symmetry mentioned in the introduc- 
tion acts on the double coset space by reversing the sign 
of X\. 

We next consider M as a product of m = L/a elements 
of G: 

Mm = T m T m -i . . .T\. (13) 

This discrete form could represent a discrete ID model 
with lattice spacing a, or could be used to take a con- 
tinuum limit a — > with L fixed. In order to increase 
the length from m to m + 1, we multiply M m on the left 
by T m+ \. In general, we take Tj to be a random vari- 
able, in which case we obtain a stochastic process. If we 
assume that TjS at different j are statistically indepen- 
dent, and that all Tj are identically distributed, then we 
obtain a process whose dynamics is "time" -independent, 
where by "time" we mean position j along the wire (this 
terminology will be used consistently throughout this pa- 
per). Also, the time-evolution process is invariant under 
right multiplication by any element of the full group G, 
in the sense that (T m+1 M)g = T m+1 (Mg) for g £ G, 
and in particular under right multiplication by any ele- 
ment of K. Thus it can be viewed as a stochastic pro- 
cess on the coset space G/K (and this is the point of 
view most often taken in the literature, see 

Refs.EllH- 

But it is not in general invariant under left multiplica- 
tion even by elements of K, unless the distribution of Tj 
is invariant under Tj — > kTjk~ l for all k 6 K and all j. 
Unless this assumption is made, we do not obtain a well- 
defined stationary stochastic process on the double coset 
space K\G/K, because the effect of left-multiplication 
by T m+ i on a double coset KgK (where each K stands 
for all members of K) depends on which element of the 
K is chosen at the left. 



It is generally believed, however, that even in a generic 
situation, due to the compactness of the subgroup K, 
the diffusion in the angular variables will lead eventu- 
ally (for long wires) to a uniform distribution on K 
(an ergodicity assumption). Therefore, one can use K- 
invariant distributions for Tj to analyze localization prop- 
erties in long wires, without loss of generality. Together 
with the if-invariant initial condition M = I for zero 
length, we can view the process as taking place on the 
double coset K\G/K. All information necessary to ob- 
tain the probability distribution for Q is contained in 
the joint distribution function of the radial coordinates 
P{x\, ...Xn;L) = P(x;L). It depends on the length of 
the wire L, and our goal is to find the equation govern- 
ing the evolution of P(x; L) when L increases. 

Each T m+ i is the transfer matrix of a thin slice, and as 
such is best parameterized as T m+ \ = exp A. Here A be- 
longs to the Lie algebra of G, Lie(G), and the probability 
distribution for A is assumed to be concentrated near zero 
since the slice is thin. Also, due to the invariance of the 
distribution, A may be assumed to be in the subspace V 
of Lie(G) which appears in the infinitesimal version of the 
Cartan decomposition Lie (G) = Lie(K)+V . (An element 
of V can be diagonalized using elements of Lie (if), as we 
used earlier, but we do not do this here.) Let us illustrate 
this point for class BD. In this case Lie(G) = so(N,N), 
Lie(-ftT) = so(iV) + so(N), and the subspace V consists of 
matrices of the form 

^=( e °o). <») 

where 9 is a real N x N matrix. This form for A corre- 
sponds to motion in G/K starting from the point M m K, 
and relative to such an origin the matrices 9 parametrize 
G/K. For class Dili, A = 9 is a purely imaginary anti- 
symmetric 2N x 2N matrix. 

It will be useful to examine the case N = 1 for class BD 
explicitly as motivation for the symmetry arguments that 
follow, and because it plays a central role in the analysis 
even for general N. In this case, G = SOo(l,l), K is 
trivial, and M is described by a single real number X\, 
which is easily seen to be X\ = ^ ■ 9j . Thus the evolution 
process is a random walk on the real line parametrized 
by x\ or 9. Clearly, if this random walk has net drift of 
either sign (i.e. the mean of 9j is nonzero), then X\ goes 
to plus or minus infinity with high probability, and the 
conductance becomes exponentially small, indicating lo- 
calization. If the mean of 9j is zero, the distribution of x\ 
will have width of order \fh (we assume that the distribu- 
tion for 9j has finite variance). In this case, the behavior 
is of the same form as the critical behavior for class BD, 
mentioned in the introduction 4 . This behavior will be 
considered in more detail below; here we wish only to 
draw attention to the central feature, which is that there 
is an obvious symmetry operation on SOo(l, 1), given by 
reversing the sign of x\, and that the drift breaks this 
symmetry. This symmetry is a discrete operation, and 
cannot be obtained by using the left and right action of 
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SOo(l, 1) on itself. There is a similar analysis for N = 1 
in class DHL SO(2,C) can be represented by the group 
of nonzero complex numbers z under multiplication (in 
fact, the eigenvalues of M are z, z^ 1 ), K = SO(2) is rep- 
resented by the complex numbers of modulus 1, and the 
additional symmetry is given by z — > z . Localization 
occurs if z drifts to either or oo; here x\ = In \ z\. 

As an aside, we see from this analysis that the dou- 
ble coset spaces for classes BD and Dili for N = 1 are 
the same. In fact, the symmetric spaces G/K for the 
two cases are also the same space: both are simply the 
real line, with translation symmetry along the line (in 
suitable variables). For Dili this relies on the fact that 
SO(2,C) is Abelian, and so K\G/K and G/K are the 
same. Moreover, there is a similar symmetry discussion 
for all three chiral symmetry classes (orthogonal, unitary, 
and symplectic, or classes BDI, AIII, and CII respectively 
in the AZ classification), and for these the double coset 
and symmetric spaces are also the same space (the real 
line) for N = 1. Therefore, at least for properties that 
can be expressed in terms of the transfer matrix, the 
critical behavior that occurs in each of these cases is the 
same. While this observation may seem trivial at this 
stage, it will become much more significant once we see 
that the N = 1 cases control the long-distance behavior 
of the critical systems at N > 1 also. Such equivalences, 
which rely on well-known isomorphisms of Lie groups of 
small rank, also occur in some other cases of the trans- 
fer matrix approach to disordered wires. In Section IVII 
below, we establish mappings between the N = 1 cases 
in the above five classes directly, in terms of hopping 
Hamiltonians. 

It is natural to ask whether there are any analogous 
additional symmetries for classes BD and Dili in the 
cases with N > 1. It turns out that there are. In 
both cases, there is a larger group that contains G as 
a connected subgroup: these are 0(N, N) and 0(2N, C). 
0(N, N) has four connected components, including two 
in which the elements have determinant —1, and the 
part SO(A, N) with determinant +1 has two connected 
components, one of which is SOo(iV, N). 0(2N, C) has 
just two components, with determinants ±1 respectively. 
Certain subgroups of these act on our groups G, map- 
ping it into itself. In particular, as we wish to leave the 
identity clement fixed, there is the action by conjugation, 
g — ► hghr x , where g £ G, and h is in the larger group. 
Examples for h of the desired form can in fact be found 
by choosing h in 0(A) x 0(A) or 0(2N), for class BD or 
Dili respectively; in both cases, such an h has the desired 
effect if it has negative determinant. Then the product 
of any two such h lies in the subgroup S[0(A)xO(A)] 
or SO (27V), respectively. (In the special case of class BD 
with N odd, we can take h to be A, which commutes 
with any element of K.) For the action on the double 
coset space, we can compare this discussion with that of 
the Weyl group, and we see that we now have the set of 
all permutations and sign reversals of the Xi, so that in 
terms of the Weyl chamber G, the additional symmetry 



is simply to reverse the sign of x\. Mathematically, what 
we have described is an outer automorphism (or involu- 
tion) of the reduced root system of G/K (outer means it 
is not connected to the identity). As these are of type Dn 
in both cases, and the outer automorphisms are known to 
be related to the existence of symmetries of the Dynkin 
diagram, we do find such a symmetry. [There are similar 
additional symmetries for each of the chiral symmetry 
classes. In these cases, the operation that sends Xi to 
— Xi for all i is such a symmetry, an outer automorphism 
of the reduced root systems, which are all of type An-i 
for these classes. However, these symmetries cannot be 
written as conjugation by elements of a covering group 
of G, unlike the BdG cases above.] 

We saw in the examples with TV = 1 that if the distri- 
bution of the random transfer matrices Tj for each slice 
is invariant under the additional symmetry (as well as 
under conjugation by K), then localization would be pre- 
vented. As localization requires that all Xi including x\ 
go off to ±oo, it is now clear that the same will be true 
for TV > 1. On the other hand, if the distribution is 
perturbed so as to break the symmetry, then localization 
may set in. [For the chiral classes, there is a similar effect: 
if the symmetry under xi — > — xi for all i is preserved, 
then localization does not occur (the system is critical) 
if and only if N is oddSS.] We must now determine the 
form of probability distributions that satisfy the symme- 
try requirements. 

A popular, and usually effective, choice is to take the 
distribution for the matrix 9 in each time slice to be Gaus- 
sian, 

P[A] ocexp ^-^MtrA 2 ^ , (15) 

where G2 is a constant that can be chosen to simplify 
some later equations, and t is the "mean free path" . For 
a Gaussian, only the second cumulant is nonzero. For 
class BD, it takes the form (for N 7^ 2) 

[&a 1 b 1 &a 2 b 2 ] c = c 2^ai a 2 ^bi b 2 ; (16) 

where C2 = a/4G2^. For a generic AT-invariant distribu- 
tion, higher cumulants are also nonzero, and should be 
expressed in terms of invariant tensors of K. Let us il- 
lustrate this for class BD. In this case, for an element k 
oiK = SO (TV) x SO (TV) we write 

*-(*£)• (17) 

where ki £ SO(iV). Then as K acts on T m by conjuga- 
tion, Tj — > kTjk^ 1 , this gives 

e -> Mfcr 1 = fciflfcj, (is) 

and we see that the second cumulant 1|16|) is expressed in 
terms of SO(A)-invariant Kronecker delta in the left and 
right indices of 9, so is invariant under SO (A) x SO (A). 
Moreover, it is invariant under 0(A) x 0(A). Invariant 
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tensors under 0(AT)xO(AT), acting in this fashion on 9, 
all have a similar form involving Kronecker 6s in the left 
and right indices. (The cumulants 



\9a 1 b 1 @a 2 b 2 ■ ■ ■ @a 2k b 2 



(19) 



must also be invariant under permutations of pairs a^, 
bi.) The number of such invariant expressions increases 
rapidly with the degree of the cumulant, and the 2fcth 
cumulant contains a distinct parameter for each one. 

However, there is also another invariant tensor of 
SO (AT), which is not invariant under O(N); the 
completely antisymmetric Levi-Civita tensor £ ai03 ,.. ew . 
Hence, there can be a term in the iVth cumulant of the 
form 



c M^a 1 a 2 ...a N Sb 1 b 2 ...b N • 



(20) 



for class BD, and 



P[A] oc exp 



C 2 £ 



trA 2 



C' 2 t 



^ ^ £a 1 a2a 3 a4^a 1 a 2 ^a 3 a 4 



ctx ,...,0.4 



(23) 



for class DHL Thus also for N — 2, one might expect 
that localization can be induced by introducing the dis- 
tributions with nonzero C' 2 (the corresponding effect on 
the DOS was observed in Ref. 0). 

In the rest of this paper, we analyze explicitly the role 
of symmetry breaking by the probability distribution for 
Tj, within the FP approach, for all N, which means that 
in general the FP equation has to be generalized to incor- 
porate non-Gaussian effects. We also perform an analysis 
(for class BD) within the supersymmetry approach. 



For even N, the iVth cumulant will contain contributions 
of both types. At higher degrees, there are also invari- 
ants that are products of both one or more Levi-Civita 
tensors and Kronecker 6s. Under the additional symme- 
try, which corresponds to acting on 9 with k\ and k 2 , 
where det k\ det k 2 = — 1, the Kronecker 6 function is in- 
variant, but the Levi-Civita tensor reverses sign. Thus 
the cumulant (|2*U|l is the one of lowest degree that (if c' N 
is nonzero) breaks the additional symmetry, but not that 
under K . 

For class Dili, the matrix 9 is imaginary and anti- 
symmetric. K = SO(2N) acts on it by conjugation. 
Consequently, there are A"-invariant cumulants contain- 
ing Kronecker 6s in arbitrary pairs of indices (the lowest 
being [9 aia2 9 a3a4 ] c = -c 2 (6 aia3 5 a2ai - 6 aiai 5 a2a3 ) with 
c 2 = a/4C 2 £), which are invariant under 0(2N), and the 
lowest-order form that is invariant under SO(2N) but not 
under 0(2N) is now 



CatE 



N^aia 2 ...a 2N i 



(21) 



which is again Nth order in 9. 

Thus in both classes BD and Dili, for A^ > 2 the use 
of a Gaussian distribution loses the symmetry-breaking 
aspect of the generic distribution, and introduces an ad- 
ditional symmetry into the evolution of the probability 
distribution with the length of the wire. On the other 
hand, for A" = 1, the cumulant that breaks the addi- 
tional symmetry is degree one, and corresponds to the 
nonzero mean or drift for 9. For A~ = 2, the symmetry- 
breaking term is degree two, the same as that obtained 
from the invariant Gaussian. Both terms in the cumu- 
lant can be obtained by using a more general form of 
Gaussian distribution^, of the forms 



P[A] oc exp 



C 2 l 



txA 1 



CU 



E 

ai,a,2 
b l i b 2 



(22) 



III. GENERALIZED FOKKER-PLANCK 
EQUATION 

In this section we first briefly review the main aspects 
of the continuum limit or FP approach for the BdG sym- 
metry classes. Then for classes BD and Dili we argue 
that while, for well-behaved distributions, the cumulants 
that break the additional symmetry are lost in the usual 
continuum limit, in general for the analysis of the long 
time behavior they should be included as higher-order 
differential operators within a generalized FP equation. 
Then we show that at long times, where the problem re- 
duces effectively to the AT = 1 case, these terms reduce 
to the drift term and produce localization. 

If one uses a Gaussian distribution of the form i|15|) 
independently for each Tj, then in the continuum limit 
a — > one obtains a heat (diffusion) or FP equation 
on G, which can be reduced to one on the coset space 
G/K (the derivation is reviewed below for a more general 
distribution of Tj). In view of the invariance under the 
left action of K of both the stochastic process on G/K 
and the initial condition, the probability distribution on 
G/K depends only on the radial coordinates a;,-, and the 
process can be viewed as one on the double coset space 
K\G/K for the joint distribution P(x; L) of the Xi (using 
the initial condition P(x; 0) = IX ^(^i))- For the BdG 
classes, this FP equation takes the form (first studied in 
the physical context in Ref . 0) : 



JV 



dL 2j£ ^ dx z dx t 



(24) 



N N 

J(x) = J[ sinh ro " 2 Xl JJ JJ sinh m °(^ ± a*), (25) 

i—l i<j ± 

where mj and m are non-negative integers (multiplicities 
of reduced roots) which depend on the symmetry of the 
problem (m = 1, m/ = for class BD, and m = 2, 
mi = for class Dili), and 7 = m (N — 1) + mi + 1. 
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[To obtain this form, we have now fixed the previously- 
introduced parameter C 2 in terms of 7, by 7 = AC 2-] The 
differential operator 

i— 1 

which appears in the right-hand side of Eq. (|24|) is re- 
lated to the part A2 of the Laplace-Beltrami operator on 
the space G/K that is independent of the angular co- 
ordinates (for the information on analysis on symmetric 
spaces, see Refs. I21U24I: relevant facts are also given in 
Appendix EJ : 

D 2 = JA 2 J~ 1 . (27) 

The factor of J arises because P is the probability den- 
sity on K\G/K, not on G/K. It is the Jacobian that ap- 
pears in the G-invariant integration measure over G/K 
in terms of the radial coordinates Xi and angular coor- 
dinates in K . For classes BD and Dili, J and D 2 are 
invariant under x\ — > —x±. 

Pertinent aspects of the solution of the Fokker-Planck 
equation are reviewed in Appendix [Bj In particular, the 
solution has a simple form for long wires when L 3> ~f£: 

P(x;L)n ^jLy^flexp^-^-m) 2 ), (28) 
i=l 

£i = ll{m (i-l) + m i y 1 . (29) 

The asymptotic form l|28|) is a product of Gaussian fac- 
tors, one for each radial coordinate. This implies that 
for large values of L/jl, the radial coordinates Xi are 
with high probability well separated compared with their 
fluctuations. This behavior of Xi is stated precisely in 
the Oseledec theorem^, which says that the xi become 
statistically-independent and Gaussian-distributed with 
means and variances proportional to L: 

(x^ = L/ii, var x { = L/j£. (30) 

In these relations the quantities £j are called localization 
lengths, and their inverses are the Lyapunov exponents. 

Here we note that the easiest way to understand Eqs. 
(|28H3U[> is to take the limit of large and widely sep- 
arated Xi directly in the Fokker-Planck equation. As 
explained in Appendix El in this limit the expression 
for the Jacobian J(x) simplifies, so that d\nJ/dxi — > 
2(m {i — 1) + to/) , leading to the asymptotic form of the 
Fokker-Planck equation: 

dP 1 A fd 2 P . . \ dP\ , s 

M=w^\^r 2{mo[l - 1)+mi) ^j- (31) 

This is a simple diffusion equation with a constant dif- 
fusion coefficient and a drift, and its solution (using the 
same initial condition as before) is Eq. (|28|) . 



The Landauer conductance JHJl of a long wire is dom- 
inated asymptotically by the smallest radial coordinate 

Xf. 

g~d/cosh 2 Xl . (32) 

Thus, the leading behavior of both the mean and the 
typical values of Q is determined by the longest of the 
localization lengths £ = £1 = J^/mi, corresponding to the 
smallest Lyapunov exponent. In most symmetry classes 
£1 < 00, and the leading behavior of Q is exponentially 
small, and this is referred to as strong or exponential 
localization. 

A distinctive feature of the FP equation for classes 
BD and Dili, first investigated in the physical context 
in Ref. 0, is that the root systems of the correspond- 
ing symmetric spaces do not have long roots: mj = 0. 
Consequently, there is no drift term for x\ in Eq. I|31(l . 
which results in x\ performing a random walk about zero. 
(We now understand this behavior as a consequence of 
the additional symmetry of these two double coset spaces 
and of the FP equations (|24|l under x\ — > — zi.) Then 
the smallest Lyapunov exponent l/£i vanishes, leading 
to the absence of exponential localization in long wires. 
In this case, we should be a little more careful about the 
treatment of x\ as L becomes large. We wish to assume 
that all \xi\ are well-separated, so that x%, . . . , xn are 
large, but not necessarily that \x\ \ is large. For the cases 
of interest, with to/ = 0, we find that J is asymptotically 
independent of x\ , and that Eq. I|31|l is still valid in this 
regime (even at small x\). The behavior of the conduc- 
tance Q in these cases is characterized by a very broad 
distribution, with an algebraic decay of the mean and the 
variance, while the variance of In £7 is of order IX: 

V ttL S 

<'" e >~-V5' varln s ~ 1 ^' (33) 

Note that the mean conductance in this situation is 
"super-Ohmic" , since Ohmic conductance decays as L^ 1 . 
The wire with this behavior of the conductance was called 
critical in Ref. 0. 

The preceding FP analysis was for the a — * limit of 
a Gaussian distribution invariant not only under K but 
under the additional symmetry. For general N , especially 
N > 2, we wish to include the effect of a non-Gaussian 
distribution for the Tj into the analysis, since for N > 2 
these are the effects that break the larger symmetry down 
to K and may lead to localization. (The following analy- 
sis also includes the cases N = 1, 2, in which these effects 
occur already at the level of general A"-invariant Gaus- 
sian distributions.) If we wish to cast this into the very 
convenient form of a differential equation in continuous 
time, then we must first consider the a — > limit for the 
more general situation of an arbitrary distribution for 
each time step. (This may be a suitable place at which 
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to mention that the use of discrete but statistically in- 
dependent increments Tj can be used to model not only 
processes that are actually discrete, but also those that 
are continuous but not (S-correlated in time, using a of or- 
der of the correlation time, provided the correlations are 
not long range. In practice, the correlation time might 
be of about the same order as the mean free path. As 
we will see, it is in these situations that the symmetry 
breaking effects for N > 2 can show up.) 

In the standard derivation of a FP equation, there 
is some distribution for each (independent, identically- 
distributed) increment of the random variables being con- 
sidered. In the following, we make the assumption that 
this distribution has finite moments of all degrees. The 
change in the probability density for the variables at time 
m + 1 from that at time m at any point is given by the 
master equation, and contains two terms, one for reduc- 
tion of the probability by transitions to all other values, 
the other for the increase by transitions from other val- 
ues (we give only a sketch of this derivation, as it can be 
found in many textbooks, and adapted to the present case 
without difficulty). Since the distribution of increments 
is narrow, the increments in the variables are typically 
small, and in the second term the probability distribution 
for time m can be Taylor expanded in the increments, 
about the values at time m + 1. The convolution inte- 
gral over these increments then reduces in each term to 
a multiple partial derivative of the product of a moment 
of the distribution for the increment and the probabil- 
ity distribution. The zeroth-order such term cancels the 
first term in the master equation. The result, in which 
we have not yet taken any limit a — > 0, is a generalization 
of the FP equation which is a first-order finite difference 
equation in the time variable m, and all orders in partial 
derivatives in the variables. Such an expansion for a fi- 
nite increment in the variables in terms of derivatives of 
the distribution is known as a Kramers-Moyal expansion, 
see, for example, Ref. |2(J 

In the present case, using a A-invariant distribution of 
the increments (described by the matrix A) , each term is 
a differential operator on G/K that is invariant under the 
left action of the whole of G. Each operator of order k has 
a coefficient that is related to a kth moment of the distri- 
bution of increments. By approximating the time step of 
length a by a first derivative, we obtain a generalized form 
of FP equation for the probability distribution P(9; L) on 
G/K, where the variables 9 introduced earlier also serve 
as the coordinates on G/K , with 9 = being the ori- 
gin, which corresponds to the transfer matrix M being 
the identity in G. The G-invariant differential operators 
on symmetric spaces have been extensively investigated, 
see for example Ref. Because sums, differences, and 
products of invariant operators are also invariant, the G- 
invariant operators form an algebra, which turns out to 
be commutative, and in fact all the operators can be ob- 
tained as linear combinations of products of a finite set of 
operators. For both classes BD and Dili, these operators 
can be labelled A2, A4, . . . , A2at_2, each of which has 



order given by its subscript, and A' PN , which is of order 
N. Thus the generalized FP equation can be written as 

ni,...,n N ' 

(34) 

where the coefficients m raii ... jraw are essentially moments 

of degree fc to t = J2k=i + Nun of the increments 9, 
except for mo,...,o which is zero due to probability con- 
servation. [Strictly, these are the moments of 9 minus 
its mean, except for the coefficient of the first-order op- 
erator, if any, which is the mean of 9. But the mean is 
zero except in the N = 1 case, which as we saw can be 
handled directly.] The sum is over nonnegative n^. 

The form of the differential operators is easy to under- 
stand if we examine the region near the origin. For class 
BD, invariant differential operators can be constructed 
from d/d9 ab . Invariance under K = SO{N)x SO(iV) 
implies that near 9 = 0, the above operators reduce to 

A — d d d d 

ajT^a^atb! d9 ai b 2 99 a2 b 2 d9 ak b 1 ' 

b 1 ',...',b k 

~ £ a 1 ...a N £b 1 ...b N -£- ""ah ' (^6) 

a x ~/o. N 09 ai b 1 09 aN b N 

i>i i>jv 

Here we see that the first N— 1 operators resemble traces, 
and were constructed using Kronecker 5s, while the last 
used the Levi-Civita tensor e. Hence, the first N — 1 
are invariant under the larger symmetry [the action of 
any element in 0(N)x O(N)}, while A' pN is odd under 
elements of 0(N)x O(N) that have determinant — 1. For 
class Dili, one has, more simply, 

A' - l — V s — 9 f38l 

Oi a 2 N uu aia,2 uu a 2 N-ia2N 

with similar transformation properties, this time under 
0(2 AT). We emphasize that in both cases these forms are 
only valid for small 9, and that for general positions on 
G/K the expressions become more complicated, in such a 
way that the operators are invariant under G (although, 
in local coordinates around any point, they would take 
the same above form, due to the symmetry of the sym- 
metric space). 

The operator A2 is the Laplace-Beltrami operator on 
G/K in both cases. For the special case N = 2, A' pN 
is also second order. It is a theorem that for an ir- 
reducible symmetric space G/K, the Laplace-Beltrami 
operator is the unique second-order invariant differen- 
tial operator, up to multiplication by a constant. For 
N = 2, the symmetric spaces for transfer matrices in 
classes BD and Dili cease to be irreducible. This is ba- 
sically because of the isomorphisms of Lie algebras so(4) 
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= sl(2)xsl(2), and the symmetric spaces likewise be- 
come direct products^. There is then a distinct Laplace- 
Beltrami operator for each factor in this product (each 
operator involving only derivatives along the directions 
lying within the factor space). The same essential point 
was noticed earlier in Ref. 0- [Incidentally, the transfer 
matrix spaces for all three chiral symmetry classes for 
N > 1 are also reducible, and contain the real line as 
a direct factor (the position on this line is described by 
Y^i x i)- Consequently, for all three cases, there are two 
invariant second-order operators (as pointed out for the 
orthogonal and unitary cases in Ref. l27|) . and a drift term 
along the real line (i.e. a first-order operator) is also in- 
variant under the simply connected group G, but breaks 
the discrete symmetry which sends 9 — ► —6, or in terms 
of the radial variables, Xi — > —Xi for all i. In a chain with 
off-diagonal randomness in the nearest-neighbor hopping, 
the drift term is produced by staggering the magnitude 
of the mean hopping 20 ^.] 

Due to the invariance of the evolution and initial con- 
ditions under the subgroup K, only the radial parts of 
the operators A 2 k and A' pN will be needed in the fol- 
lowing. In terms of the radial parts as defined in Ap- 
pendix [X] A' pN = N\A PN . The angular dependence of 
P(0; L) is trivial, and we can reduce the equation to 
one for the probability density in only the Xi variables, 
P(x; L). Because P is a probability density, the relation 
is P(x;L) = JP(9;L). The generalized FP equation for 
P(x;L) takes the same form as Eq. (I34f) . but with D 2 k 
and N\D PN in place of A 2 k and A pN , where 

D 2 k — JA 2 kJ 1 , 

D PN =JA PN J-\ (39) 

We can now describe the fate of the terms of more than 
second order in the generalized FP equation, as a — ► 0. 
For any probability distribution of Tj which has finite 
second cumulant with coefficient c 2 for 9, we define the 
mean free path £ through Eq. , where C 2 depends on 
the symmetry class, and is chosen so that the coefficient 
of the second order operator in the FP equation takes 
the standard form, Eq. We are only interested in 

a < £. If has nonzero mean, which can occur only for 
N = 1, then we take this mean to be a number times 
a/ ("/£). The moments m ni! ... )njv (other than the mean) 
with total degree fctot = ^2k=i ^krik + Nun are then 
of order (a/j£) ktot ^ 2 . If we take the limit a — ► with 
7^ fixed, then all operators of higher order than second 
drop out, and the usual FP equation is obtained. For 
TV > 2, all terms that break the additional symmetry 
are lost in this limit. For N = 1 and 2, a drift term, 
and an additional second-order term respectively survive 
in the limit, and these do break the symmetry unless 
their coefficients are tuned to zero. It has been noticed 
previously that localization ensues when the drift term 
is present for TV = 1, as discussed abovei^, and that the 
non-critical form of the DOS, Eq. JIJ, results when the 
additional second-order term is included 9 for N = 2 . 



[Similarly, for the three chiral classes, the N transitions, 
at each of which the behavior is critical, can be reached 
within the FP approach by tuning the coefficient of the 
drift termSSl.] 

It is instructive to realize that the vanishing of all 
higher-than-second-order terms in the limit a — > is con- 
nected with the central limit theorem, and also to inter- 
pret it in terms of the renormalization group (RG). When 
the typical values of 9 a b for each time step are small (com- 
pared with 1/y/N), then in the product of Tj over any 
time interval that is not too long (less than about I — not 
j£, as we explain below), to a good approximation the 
9s simply add (as matrices). G/K is a smooth manifold, 
and has a finite radius of curvature at any point. Over 
this time interval, the random position on G/K has been 
displaced by an amount less than the radius of curva- 
ture. Hence we are essentially looking at a random walk 
in a flat space. The central limit theorem for sums of in- 
dependent, identically-distributed random variables then 
tells us that the sum becomes Gaussian distributed in the 
limit of many variables, provided the mean and variance 
for each variable are finite. This means that the proba- 
bility density for the sum effectively obeys a FP equation 
with no higher order terms over a long enough time in- 
terval. In this flat space case, this can be understood in 
terms of the generalized FP equation as saying that as the 
time gets longer, the distribution spreads, and the values 
of the higher derivatives become smaller relative to the 
second derivative, and so can be neglected. In renormal- 
ization group language, the terms of higher order than 
second are irrelevant (in this flat space regime), as they 
become less important at larger time scales. However, if 
we reverse this process and go to shorter time scales, the 
higher derivatives get larger. If we take the full general- 
ized FP equation literally and use it for the evolution of 
an initial S function for arbitrarily short times, then it 
ceases to make sense. Of course, the derivation from dis- 
crete time steps implies that the equation for time step 
a certainly cannot be used for times less than a. It effec- 
tively describes the process "at the time scale a". The 
fact that the irrelevant operators have coefficients that 
contain positive powers of a (the short-time or ultravio- 
let cutoff) is typical for irrelevant operators. As always 
in the RG, "relevant" and "irrelevant" refer to linearized 
scaling behavior near a particular fixed point. In the 
present case, the fixed point is the ballistic fixed-point, 
which is described by the FP equation, or by using a 
Gaussian distribution for the increments and taking the 
a — > limit. The a — > limit can be taken only if we are 
prepared to drop the effects of the irrelevant operators. 

We now examine the behavior at longer times. As 
its name implies, the ballistic fixed-point describes very 
short times where there is no scattering. The mean free 
path £ is the time scale at which scattering, which cor- 
responds to curvature of G/K, first shows up. We are 
interested in a small, a<f, and all the coefficients of the 
irrelevant operators can be thought of as small, though 
we do not need to set them to zero. Over these time 
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intervals, the central limit theorem applies, and the FP 
equation gives a reasonable description of the evolution 
ofP(a;;L). For L < £, all small, and Q ~ dN . We 

should note that the squares of the XiS are the eigenval- 
ues of 98 T , and so have typical values of order N times 
those of the squares of elements Q a b- As the time L passes 
through the localization lengths £n < 6v-i < •• • < £2, 
typical values of the corresponding Xi pass through the 
vicinity of 1 and begin to grow linearly (the smallest such 
length £jv is close to £). Their contribution to Q there- 
fore drops, and we are in the Ohmic regime in which the 
mean and typical Q are of order dN£/L (Nl/L is the 
number of channels not yet localized at time L). Here 
there are also fluctuations in Q. This crossover behav- 
ior is "universal" , even at finite N, in the sense that 
the irrelevant operators can be dropped; the universal- 
ity arises, as usual in the RG, from the RG flow passing 
close to a fixed-point, in this case the ballistic fixed point 
which is reached as L — > provided the irrelevant oper- 
ators are dropped. However, this universality may not 
be of great interest physically. If one also takes the limit 
N — > 00 with 7^ fixed, so I — > 0, before taking L — > 0, 
then one arrives at a "diffusive" fixed point. That is, if 
N is large but finite, then the RG flow from the ballistic 
fixed point approaches close to the diffusive fixed point. 
This universality in the Ohmic regime when the diffusive 
fixed point is approached as L — > is what is usually 
meant when the term "universal" is used in mesoscopic 
physics for the behavior of the conductance in a ID sys- 
tem. The diffusive fixed point will not concern us here, 
but we point out that the higher-order differential oper- 
ators are also irrelevant at this fixed point, and so are 
dropped to obtain this universal behavior also. Then not 
only the Ohmic regime, but the whole crossover that be- 
gins from the diffusive fixed point at L — > is universal, 
and can be studied as in Ref. 0. 

When the time L passes approximately the inverse of 
the last nonzero Lyapunov exponent, that is L ~ £2 — 
"f£/m , then only Xi is still not growing linearly. If one 
uses the FP Eq. (|24|l . then this is where the asymptotic 
behavior termed critical sets in. Thus the Ohmic regime 
is a crossover from the diffusive fixed point to a critical 
fixed point, and £2 is our estimate for the crossover time, 
£x ~ £2 ~ j£. In most other symmetry classes, local- 
ization would set in at times larger than £1, so that the 
localization length £ is about in agreement with RG 
analysis in ID. In the present case, the analysis of the RG 
is similar, but the crossover is to a critical rather than a 
localized fixed point. 

We are concerned with the effect of the symmetry- 
breaking operators in the generalized FP equation on the 
critical fixed point. We assume that the coefficient of the 
leading such term (D PN ) is small, but not zero. Then it, 
and all other terms of orders higher than second, can be 
dropped when studying the crossovers from ballistic to 
diffusive to critical. The next question is whether they 
are important (relevant) at the critical fixed point. 

The main points can be understood by retaining in 



the right-hand side of the generalized FP only the terms 
containing the differential operators Di and D PN , so that 



dP{x;L) 
dl 



(40) 



1) c' N ). Because we are in the 



(note that TOc^o,...,! = 
regime where all Xi except x\ are large, and because the 
asymptotic behavior for c' N = is dominated by xi, we 
will attempt to "integrate out" these X{. In Appendix IA1 
we determine explicitly the asymptotic forms D2k and 
D PN of the operators D 2 k and D PN for classes BD and 
Dili, in the regime \x±\ <C X2 <C ■•■ -C xn (but no 
other assumptions on x\). Then we define the reduced 
probability density Pi(x\,L) for x\ by integrating Xi, 

1 2 over the range \xi | ^ X2 ^ ■ ■ ■ ^ £jv- An equation 
for Pi can be obtained by integrating the generalized FP 
equation over the same range. By using the asymptotic 
form of P(x;L), Eq. (|28[l . (in the derivation of which, 
the effects of c' N and other higher-order operators are 
negligible), we find that in the D's, the d/dxiS for i ^ 

2 can be replaced by zero to obtain the corresponding 
differential operators in x± only, which act on Pi. In 
particular, D PN becomes oc d/dxi, and still breaks the 
symmetry x\ — > —x\. The resulting equation (for either 
class BD or Dili) is 



dPi 
dL 



1 d 2 Pi dPi 



2£ x dx{ dx\ 



Here we used 



Cx = it, 



N- 



\N - 1)! 



-N 



(41) 

(42) 
(43) 



Equation l|41() has no explicit ./V-dependence, and in par- 
ticular is precisely the form obtained for N = 1, on 
dropping operators higher than second order. Also, this 
justifies the use of the term "critical fixed point" when 
k = 0, since the equation allows a scaling analysis, with 
x\ ~ (L/^x ) 1 / 2 . This fixed point reproduces the critical 
asymptotic behavior of Q or hit/, of course. 

Hence, the symmetry-breaking terms produced by dis- 
order that is not invariant under the additional symme- 
try are relevant at the critical fixed point, and there is 
a phase transition at k — between two phases. In the 
RG, such operators that are irrelevant at a fixed point 
in the ultraviolet (such as the ballistic or diffusive fixed 
points), but relevant at a fixed point at a larger time 
scale (in the infrared) are called dangerously irrelevant, 
because even though they are irrelevant at the ultravio- 
let fixed point, they cannot be set to zero without some 
important physics being lost. For k 7^ 0, Eq. (|41|> shows 
that the Lyapunov exponent for x± is so that 



(44) 
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describes the localization length when the typical (i.e. 
high-probability) behavior of — In Q is growing linearly 
as 2L/£. The behavior at the transition is universal 
and described by Eq. (|41|l provided L and £ ^> £ x , 
even if N is finite or small so that the crossover from 
the diffusive to the critical fixed point is not fully uni- 
versal. In this regime, the probability density of y = 
(— i ln0)/(L/^ x y/ 2 approaches a universal function 
of y and £ x L/£ 2 , because lnC* — ln(4d) is approximately 
— 2|xi|, except near Q ~ d, but that region has negligible 
probability for large L. The precise result is 



P(y) 



exp 



1/2' 



exp 



1/2 



(45) 



The regime of linear growth of the typical value of — In Q 
sets in at the time scale L ~ £ 2 /£x , which is much larger 
than £ under our assumptions that define the scaling 
regime. At smaller time scales, it grows only as C^L/^x, 
for some constant C, as described for the critical k = 
case using the N channel FP equation. 

All the moments of Q are of order e~^ x L '/( 2 £ ) y/CxjL, 
throughout the critical regime. (This behavior differs 
from the "log-normal" behavior in a one-channel wire in 
most other symmetry classes, because usually the analogs 
of £ x and £ are proportional.) Thus, while £typ = £ = 
1/\k\ describes the rate of linear decay of the typical log 
conductance, 



point). This thinning out of degrees of freedom, a typ- 
ical method in the RG, would show explicitly how the 
RG acts in this regime. The variables Xi that are grow- 
ing linearly correspond to channels in which transmission 
is decaying exponentially, like massive degrees of freedom 
in a field theory. 

We can combine the forms obtained here for the mean 
and typical localization lengths with the values of the 
moments that would be obtained for the non-Gaussian 
probability distribution of each Oj at the lattice scale a. 
If the distribution has a fixed form that scales with the 
width ~ (a/7^) 1 / 2 , with a parameter h for the strength of 
the leading symmetry-breaking cumulant in these units 
(h stays fixed as a — ► 0), so that c' N = h(a/ 1 £) N / 2 for 
N ^ 2, then we find that 



£ty P 



1 



m^ 1 



,a 



(n — iy.\h\ 



-(JV/2-1) 



(47) 



while £mean/£x = 4(£ t yp/6x) 2 - This illustrates how the 
(dangerous) irrelevance of the symmetry-breaking per- 
turbation at the ballistic and diffusive fixed points may 
translate into very large localization lengths, especially 
as a/ji — > with N large but fixed, and £ x and h fixed. 
Note that the limit N — > oo with £ x fixed is appropriate 
for studying the diffusive regime and the crossover; the 
localization lengths appear to vanish in this limit if h and 
a/j£ are held fixed, however the analysis assumed that 
£typ ^ £x (i-e. small \h\). That is, if \h\ is sufficiently 
large, localization will occur without the RG flow ap- 
proaching the vicinity of the critical fixed point analyzed 
above. 



Cmean — 4£ /£> 



(46) 



describes the exponential decay of the mean conductance. 
Similar behavior — L/£ t y P and e _L ^ mcan can be inferred 
for the typical behavior of the logarithm, and for the 
mean behavior, respectively, of the eigenfunctions at zero 
energy, as these determine the transmission through the 
system. This behavior also applies to mean values of the 
Green's functions, as we will see in the following section. 
The mean localization length is much longer then the typ- 
ical one, because the average transmission is dominated 
by the rare events that produce \x\\ < 1. 

We can check the effect of all other terms in the gen- 
eralized FP equation in the same way. We find that in 
asymptopia, all D^k reduce to operators d 2k /dx\ k . At 
the critical fixed point, operators of order higher than 
two are again irrelevant, and the only other term that 
should be kept is a correction of second order coming 
from (D PN ) 2 , which gives a small correction to l/£ x - As 
an aside we mention that a similar analysis could be done 
in the (Ohmic) crossover regime, by integrating out only 
the subset of Xi that typically are growing exponentially 
(this means introducing a somewhat arbitrary cutoff on 
i, since there is not a clear separation of large and small 
Xi in this case, unlike in the vicinity of the critical fixed 



IV. SUPERSYMMETRY SOLUTION 

In this section, we discuss the application of a super- 
symmetry method to the problem of localization in a su- 
perconducting wire in class BD. The details of the calcu- 
lations are provided in Appendix^ and here we only give 
a summary of the results. The technique is effectively a 
finite-iV analog of the ID supersymmetric nonlinear a 
model, which arises in the N — > oo limit. The forward 
and backward directions along the wire are referred to 
here as "up" and "down" , and after averaging over disor- 
der in the subgroup K, the states are in a tensor product 
space of two "superspins" , which are the representation 
spaces of the superalgebra osp(2|2), R and R for up and 
down, respectively. After averaging, these are coupled by 
the backward scattering (9) part of the transfer matrix 
M. The usual Gaussian disorder, invariant under the ad- 
ditional symmetry, produces a coupling of the superspins 
that corresponds to the usual Heisenberg coupling famil- 
iar from magnetism. The leading symmetry-breaking cu- 
mulant produces an additional term with the coefficient 
k, given by Eq. I|43|l . When it is nonzero, it produces lo- 
calization. A dimensionless parameter |A| = £x/(2£typ) 
that enters the calculations in Appendix El is propor- 
tional to |/d. 
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The properties of the effective superspin problem are 
easily translated to the properties of the disordered wire. 
In particular, the mean localization length £ mcan in the 
wire is inversely proportional to the energy gap above the 
ground state at zero "energy" (conjugate to "time" , i.e. 
distance along the wire) in the superspin problem. The 
DOS in the wire is proportional to the expectation value 
(in the ground state) of the "staggered magnetization" 
for the two superspins, as a function of true energy. We 
have solved the superspin problem and have found its 
eigenstates and the spectrum. (We present results for N 
even only. The case N = 1 was done previously^, using 
results from Ref. |2^.) The excited states are separated 
from the ground state by a gap that turns out to be 

-^gap — 2/£ mcan , (48) 

where £ mca n is again given by Eqs. (^HlElffl} f° r m o — 
1. This implies that all moments of Q will again decay 
as ~ e _2L /^ moan (when E gap is nonzero), confirming the 
result of the previous section by a very different method. 

We also determine the low-energy DOS v(e) in the wire 
both at and near the critical point. At criticality (A = 0), 
we obtain Dyson's singularity (J2J). Close to the critical 
point (|A| <C 1) we get the divergent DOS as in Eq. JJJ 
with exponent 5 = 1 — 4|A|. At criticality but at finite e, 
the behavior can be characterized in terms of a typical 
localization length £ t yp ~ | In e| as e — > + , and a mean 
localization length £ moan ~ In 2 |e|. This distinction has 
been emphasized by D. Fisher in his study of random 
spin chains^. The same forms that we find here were 
obtained using supersymmetry techniques in Ref. 28 for 
a one-channel model in the chiral orthogonal class (BDI 
in AZ classification). We comment further on the rela- 
tion of the transitions in ID in different symmetry classes 
in Section IVll below. It is interesting that two analogous 
length scales emerged rather easily from the transfer ma- 
trix analysis (especially if one considers only N = 1). 
For the random hopping problem (class BDI with one 
channel), an approach using the transfer matrix and a 
scaling argument^ reproduced Dyson's DOS at critical- 
ity. It appears that one can think of the DOS integrated 
from to e (which can be obtained from the transfer 
matrbiS) as being of order l/£ moan by scaling, as it is 
really a mean Green's function property; then differenti- 
ating reproduces Dyson's result. 

It was emphasized by Bocquet et alJ-L that the target 
space of the nonlinear a model for classes BD and Dili is 
disconnected (it is OSp(2|2)/U(l|l) for the simplest case 
in class BD, and the bosonic sector is 0(2)/U(l), which 
has two components). Thus there can be domain walls at 
which the field jumps from one component to the other. 
In terms of superspins, the representation R splits into 
subspaces called i?o and Rn in Appendix IO which are 
irreducible representations of the superalgebra osp(2|2) 
(there is a similar decomposition of R). Under the ac- 
tion of the superalgebra, the orbits of the highest or low- 
est weights in these subspaces form supermanifolds which 
are the two connected components of OSp(2|2)/U(l|l), 



so that states in these subspaces correspond to being on 
one or other component of the supermanifold. The inter- 
esting point is that the term that breaks the additional 
symmetry produces a term in the superspin Hamiltonian 
that flips between these two subspaces, and thus acts as 
a domain wall as a function of time (the domain wall 
is simply a point in ID). Hence we can conclude that 
domain walls are essential to produce localization. This 
has been proposed befor e) 11 ' 12 , but now we have explicit 
results for arbitrary numbers of channels in ID, and in 
our framework we know the precise form of the effect the 
domain wall produces (though we have not analyzed the 
N — > oo limit). 



V. RANDOM FLUXES 

So far, our analysis has been performed in the continu- 
ous model introduced in Section \n\ in which the transfer 
matrices Tj that are multiplied together for each time 
step are in the connected group G, and the disorder is 
weak so each Tj is close to the identity. But the symmetry 
conditions that define classes BD and Dili do not impose 
this continuity. They allow the transfer matrix, or its in- 
crements Tj, to belong to a larger group G in each case, 
with G = 0(N, N) for class BD, G = 0(2N,C) for class 
DHL These groups include matrices of determinant ±1, 
and so are not connected. Such transfer matrices may 
then be expected to arise when the model is not contin- 
uous, for example lattice (or tight-binding) models. The 
negative determinant may be interpreted as saying that 
when a particle propagates around a closed loop (going 
first forwards along the wire, then backwards by a dif- 
ferent path), it may encircle a magnetic flux of half a 
quantum, often referred to as a flux of it. Thus we term 
this the inclusion of random fluxes (of it) in the model 
(it would perhaps be more accurate to call these random 
vector potentials). In this section, we analyze the effects 
of these. The analysis is straightforward in view of the 
previous results. 

It is convenient to consider the following simple model, 
that is based on the continuous one used up to now. We 
use a small density (along the wire) of random fluxes. 
That is, generically the evolution is by TjS in G with weak 
disorder, distributed as before. But at random times, we 
also insert a matrix of negative determinant. Such a ma- 
trix could be an arbitrary element of G. One may wish to 
view it as a Tj close to the identity, times some element 
with a negative determinant. As we factor out the Tj 
close to the identity anyway, we may as well choose the 
extra matrix to belong to a maximal unitary subgroup 
K of G, which contains K. Relative to the same basis we 
have used throughout, which refers to the decomposition 
of M' into in and out subspaces, the natural choice for K 
is 0{N) xO(N) or 0(2AT) for classes BD and Dili, respec- 
tively (other choices for K or K are isomorphic to these) . 
With this choice, our remarks above about picking up a 
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phase —1 on propagating around some loops describe ex- 
actly what takes place, up to the effect of multiplication 
by some element of K (if we neglect the disorder outside 
K which is small over a small time-step). 

We note that for class BD, K has four connected com- 
ponents, as may be seen by writing its elements in the 
form of Eq. (JEJ but with k\ and k 2 in 0(A). Then 
the four choices of signs for the determinants of k\ and 
k 2 label the four components. For class Dili, there are 
just two connected components of K. Because 0(A) 
is a semidirect, not a direct, product of SO(N) and Z 2 
when N is even, there is no unique (i.e. invariant under 
K) choice of negative determinant matrix in K for class 
Dili or for class BD for N even. For class BD with N 
odd, there are two possible natural choices. If we write 
elements of K in the form of Eq. I|17[l but with ki and k 2 
in 0(A), then we can take one of kx, k 2 equal to In, the 
other equal to — In- 

As time evolves (i.e. L increases), the transfer matrix 
M evolves on G by diffusion as before, but now also with 
occasional jumps between the different connected compo- 
nents of G, due to the occasional insertion of a negative 
determinant element of K into the product of transfer 
matrices for small steps. Modulo multiplication on the 
right by elements of K , M evolves on the coset space 
G/K, which has four or two connected components for 
classes BD and Dili respectively, usually by diffusion but 
with occasional jumps to another component. Modulo 
the left action of K on these symmetric spaces, evolution 
is on four (two) copies of the double coset space K\G/K 
for class BD (Dili, respectively). For class BD, the dis- 
tinction between components with the same sign for the 
determinant appears to be unimportant, so we disregard 
it. The important point is that there is a correspondence 
between points on the two copies of K\G/K, and when 
a jump occurs, it is between two corresponding points. If 
we parametrize the two connected spaces under this cor- 
respondence by the usual coordinates Xi, then the only 
difference between evolution on the two components is 
that in the generalized FP, which applies between jumps, 
the terms that break the symmetry x\ — > — x\ have op- 
posite signs on the two components. 

We may now consider the effect of the random fluxes on 
the wire over very long time (length) intervals that con- 
tain many jumps. If the jumps are correlated in pairs, 
so that the system spends most of its time on one com- 
ponent (say, the one in which M is in the identity com- 
ponent), then there will be a net effect of the symmetry 
breaking terms. In the critical regime, there will be a net 
drift term. In this case, a transition as we have described 
it in the previous sections will still occur, and we expect 
the same universal properties. This is what occurs in the 
lattice models considered in Refs. I7IIT7I 

On the other hand, if the jumps are uncorrelated, or 
more generally if equal time is spent on both components, 
then the average drift will be zero. In this case, localiza- 
tion is suppressed and the system remains critical. (We 



note that both types of behavior were observed in ID 
N = 1 models in Ref. 0.) In this case, the results can 
also be obtained by another approach. After moderate 
time intervals, the transfer matrices Tj become uniformly 
distributed over the disconnected angular group K. We 
can treat this situation as we did for the connected group 
K in the continuous model, which leads us to consider the 
coset spaces G/K and the double cosets K\G/K. These 
are connected spaces, and G/K is the same as G/K. 
The Weyl group inside K now acts on the radial coordi- 
nates Xi by permutations and arbitrary sign changes, so 
they lie in the Weyl chamber 0<a:i<a:2<---<a;Ar, 
which is the quotient of the previous C by the operation 
x\ — > —x%; this parametrizes K\G/K. In this case, the 
generalized FP equation on G/K can contain only dif- 
ferential operators invariant under K, and so A PN is not 
allowed, and its place among the generators of the alge- 
bra of invariant differential operators is taken by A2JV 
which is of order 2N. It is now clear from our preceding 
analysis that in this case there can be no relevant pertur- 
bation of the critical fixed point, and the system remains 
critical, with behavior as described in Ref. 0. 

We can also see how these effects would emerge in the 
supersymmetry framework for class BD (see Appendix lOl 
for terminology). If the angular variables must be av- 
eraged over K, then this produces a constraint that the 
states on the up (down) sites must be in R (Ro). Hence 
there can be no XH\ term in the Hamiltonian, and in 
this case all states will be extended. Thus, suppressing 
domain walls suppresses localization. (Again, this was 
discussed previously for 2D, and for N = 1 in lDiLi£.) 

For the chiral symmetry classes, only the chiral orthog- 
onal class, BDI, possesses a disconnected transfer matrix 
group G = GL(N, R), while K = 0(N), and both of these 
groups have two connected components. However, unlike 
the cases of classes BD and Dili, in this case K acting by 
conjugation on the Tj does not reverse the sign of all Xi 
(the desired additional symmetry operation). Thus the 
fluxes have apparently a trivial effect, and the discon- 
nectedness of GL(iV, R) is like the factor of 2 in the four 
connected components of 0(N 7 N) that likewise had no 
apparent physical consequences. 

VI. SUPER-UNIVERSALITY 

Our results have established universality near the criti- 
cal point, at least for the conductance and the mean den- 
sity states, in the sense that the universal properties do 
not depend on the number of channels, or the details of 
the disorder, within the symmetry classes BD and DHL 
Similar results were already known for the orthogonal, 
unitary, and symplectic chiral symmetry classes (classes 
BDI, AIII, and CII, respectively). More surprising is that 
the effective FP equation for the transfer matrix near the 
critical point, Eq. I|41|) . is the same for all five classes, 
which implies that the scaling forms for the distributions 
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both of the conductance and of its logarithm are the same 
for all five classes. Moreover, the mean density of states 
also has the same form, first found by Dyson& for one 
case, in all five cases. This suggests that there is "super- 
universality" , meaning that all five transitions are in the 
same universality class, even though they lie in distinct 
symmetry classes. That is, all corresponding universal 
properties should be the same in the critical region in 
each class. It is by now clear that in each symmetry 
class the universal properties of the transition can be ob- 
tained from the N = 1 case. It is of interest to prove the 
super-universality directly, and not only for properties 
related to the transfer matrix. Here we point out that 
this can be done within lattice models for N = 1; these 
do not include random fluxes. 

Some results in this direction exist; the equivalence of 
an N = 1 model in class BD, which possesses a transi- 
tion, with one in class BDI was pointed out in Ref. |2^. 
This can be formulated as follows within a ID tight- 
binding chain: a Hamiltonian that consists solely of ran- 
dom imaginary nearest-neighbor hopping (which lies in 
class BD with N = 1) is equivalent by a simple (disorder- 
independent) gauge transformation to a model with real 
random hopping (which lies in class BDI with N = 1). 

Now, for each of the three chiral classes, the one- 
channel Hamiltonian can be taken to consist solely of 
nearest- neighbor hopping that, when hopping to the right 
on any link, is a random positive real number (which we 
call the magnitude of the hopping) times a random el- 
ement of Z2 (i.e. ±1), U(l) (i.e. a complex number of 
modulus 1), or SU(2) (i.e. a 2 x 2 unitary matrix of 
determinant 1 — in this case only, we consider particles 
with spin 1/2), for the orthogonal, unitary, and sym- 
plectic cases, respectively. Clearly, we may perform a 
randomness-dependent gauge transformation that elim- 
inates the group-valued factors and makes the hopping 
real and positive everywhere (we assume open, not pe- 
riodic, boundary conditions). Then the spectrum of the 
Hamiltonian (for fixed disorder) is the same in each class, 
and independent of the values of the group-valued dis- 
order. The Green's functions are related by the gauge 
transformations. If we assume the magnitude of the hop- 
ping is statistically independent of the group- valued fac- 
tor on each link (and the random variables on different 
links are all independent), then gauge-invariant moments 
of Green's functions are equal for the different classes. 
These statements are independent of the distributions 
for the disorder, apart from the independence properties 
already stated. The disorder could be weak, consisting 
of small random corrections to a constant (or more gen- 
erally, staggered) real positive hopping. This model is 
useful in taking the continuum limit to establish contact 
with other models in the literature. The disorder could 
also be taken with the group-valued factors uniformly 
distributed over the relevant group, on each link. In 
this case, the "angular" (group- valued) disorder is strong, 
though the disorder in the magnitude could still be weak. 
In all cases, it follows that all universal properties at the 



transitions in these classes are the same. (It would be 
surprising if these results were not already known for the 
chiral classes.) 

This leaves only a mapping for a ID, TV = 1 model 
in symmetry class Dili into either class BD or one of 
the chiral classes to be found. It turns out that there is 
a mapping to the chiral unitary (AIII) class that paral- 
lels that of class BD to the chiral orthogonal class. We 
take a tight-binding chain with two basis states per site. 
Class Dili Hamiltonians can be considered as lying in 
class BD, but with additional restrictions. Thus we may 
begin with a purely imaginary (and Hermitian) Hamilto- 
nian on this chain, and as usual we expect that nearest- 
neighbor hopping will be sufficient for our purposes. On 
examining Altland and Zirnbaues^, or Motrunich et al£ 
[especially their equation (5)], we see (further explana- 
tion is given below) that a Hamiltonian in class Dili can 
be obtained if there are no on-site terms, and if the hop- 
ping matrix between neighbors (going to the right) takes 
the form i times a positive real number, times an orthog- 
onal matrix in 0(2) with negative determinant. By a 
disorder-independent gauge transformation, this can be 
brought to a real positive number times an SO (2) ma- 
trix. Since this is equivalent to a complex number, we 
see that the Hamiltonian is equivalent to two copies (one 
the complex conjugate of the other) of a Hamiltonian in 
the chiral unitary class. Hence our proof is complete, and 
all the universal properties of the transitions in these five 
symmetry classes are the same. 

We append here a more detailed explanation of the 
structure of the N = 1 class Dili hopping Hamiltonian. 
Let the components of the state vector be j/j, with y%j—i, 
?/2 j belonging to the same site j, j — 1, . . . , M (in 
this section, M is a positive integer). Then the Hamil- 
tonian must be an antisymmetric imaginary 2M x 2M 
matrix, i.e. an element of the Lie algebra so(2M) of 
S0(2M ). By definition-, a Hamiltonian in class Dili lies 
in so(2M)/u(M). To understand the action of u(M), we 
form the combinations Zj = U2j-i J ri]j2j- If all y,i are real 
(an assumption which does no harm), then the similar 
combinations Zj — Jfty-i — iy2j are simply the complex 
conjugates of Zj. Now a Hamiltonian that is an M x M 
Hermitian matrix acting on the vector of Zj can be viewed 
as an antisymmetric matrix acting on the y^. These 
parametrize a linear subspace of so(2M), and matrix ele- 
ments of this form are supposed not to be present. That 
means there are no onsite terms, and that any non-zero 
nearest-neighbor terms must violate the complex struc- 
ture, effectively mapping Zj to a complex number times 
~Zj+i- In the yi components, this yields the positive real 
magnitude times i times a 2 x 2 orthogonal matrix with 
determinant —1, as claimed above. The nearest- neighbor 
structure is essential to the equivalence, as the set of all 
Hamiltonians in class Dili is certainly not the same as 
that for AIII, according to the symmetry classification. 
However, there are also Hamiltonians with longer-range 
hopping that lie in both classes, and they involve the 
bipartite lattice structure. 
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We emphasize that the Hamiltonians here for N = 1 
channel in classes BD and Dili do not describe supercon- 
ductors in any simple direct way, as one hnds that TV is 
even in those models. But they are significant as we have 
shown earlier that N — 1 controls the critical properties 
in all cases. 

The equivalence of the universality classes for the tran- 
sitions is not spoiled by the fact that for classes BD and 
Dili only, one can include random fluxes that prevent the 
system from becoming localized (at e = 0) . This only af- 
fects whether or not the relevant perturbation that takes 
the system off the critical point is possible. The proper- 
ties at criticality (such as the DOS) are still the same. 

The equivalence of the N — 1 models considered here 
suggests that they may possess additional symmetries, 
especially when viewed within the supersymmetry for- 
malism. The corresponding nonlinear sigma models arc 
certainly not the same, but these correspond to the 
TV — > oo limit. Apparently, their TV = 1 analogs must be 
equivalent. One example of this phenomenon is known 29 . 



VII. CONCLUSION 

We have shown that zero-energy quasiparticle states 
in a long disordered superconducting wire with broken 
spin rotation invariance and an arbitrary number (TV) of 
channels are generically localized, both in cases with and 
in cases without time-reversal invariance. This leads to 
exponential decay of the thermal conductance of the wire 
with its length. Within a model with disorder invariant 
under the "angular" symmetry group K, there are two 
localized phases, separated by a critical point which can 
be reached by tuning a parameter characterizing the dis- 
order distribution in the wire. One of our central results 
is a universal scaling form for the probability distribu- 
tion of the logarithm of the conductance near criticality, 
Eq. H45fl . which exhibits two different length scales, which 
correspond to the localization lengths for the mean con- 
ductance and for the typical log conductance. We also 
found the mean density of states at low energies both at 
and near criticality, which agree with Eqs. CD, ©• Our 
results are in full agreement with those of Motrunich et 
al£, who considered the density of states in models with 
N = 2 or 4 channels. Our results confirm that the uni- 
versal results are independent of the number of channels. 
We showed that the universality classes for the transi- 
tions in the two symmetry classes are the same, and the 
same as those in the chiral symmetry classes. We also 
showed how the presence of random vector potentials or 
fluxes may suppress localization, leaving the system crit- 
ical. Our analysis explains why the off-critical behavior 
was not detected in some earlier work that used the FP 
eqiii. as the use of a continuum (a — > 0) limit, or of 
^-invariant Gaussian disorder, leads to the loss of the 
relevant perturbation if 7V > 2. 

One issue we have left unresolved is the nature of the 
possible multicritical point found within some class BD 



models 7 . Another more technical issue is that our analy- 
sis used the generalized FP equation, which was derived 
on the assumption that all moments of the increments are 
finite. It would be desirable to weaken this assumption as 
far as possible (as in the proof of the usual central limit 
theorem, which requires only that the second moment of 
the increments be finite). 

We want to emphasize here that our results for the 
critical regime are universal, regardless of the number 
of channels, and thus regardless of whether the diffusive 
fixed point is reached as I — > 0, which occurs only if 
N — » oo. The critical (scaling) regime is defined by the 
condition that the typical localization length £ be much 
larger than the crossover length £ x from the short length 
behavior; the critical point corresponds to £ = oo. £ x 
may be of order the mean free path £, or much larger 
~ N£ in the many-channel case in which short-length be- 
havior is diffusive (Ohmic). We used a model description 
with disorder invariant under the angular subgroup K. 
More general models do not have this property at short 
scales, but it is expected that the angular variables in the 
transfer matrix rapidly become uniformly distributed, so 
that our results should be universal for all such models, 
in the scaling regime stated. 

We also wish to emphasize that our results are generic 
for classes BD and DHL As in the three chiral symmetry 
classes (BDI, AIII, and CII), states are generally local- 
ized, but there can be phase transitions as a parame- 
ter is tuned. This critical behavior is intrinsic to these 
symmetry classes; the behavior must be described by 
a parameter that specifies the distance from the criti- 
cal point, as well as a mean free path. This may seem 
to contradict some widely-held assumptions for ID dis- 
ordered wires, that the symmetry class determines ev- 
erything, and that the presence of any other parameter 
means "non-universality" . However, in the general the- 
ory of critical phenomena, universal scaling functions of 
more than one variable are common. The necessity of 
using other parameters in addition to the conductance 
at the scale of the mean free path to describe universal 
phenomena should be familiar from the example of the 
quantum Hall effect in two dimensions, where the ad- 
ditional parameter is the Hall conductance. Symmetry 
classes are not universality classes; universality occurs 
when the RG flow passes near a fixed point, and there 
may be more than just the ballistic, diffusive, and local- 
ized fixed points within a given symmetry class. Which 
universality class occurs in a given model within a sym- 
metry class depends on other aspects of the structure of 
the model. Examples of this have been previously ob- 
served in class BD in two dimensionsi 2 * 2 ^. 
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APPENDIX A: INVARIANT DIFFERENTIAL 
OPERATORS ON G/K 

In this Appendix we gather necessary information on 
the spaces of transfer matrices for classes BD and Dili, 
viewed as symmetric spaces G/K. A standard detailed 
reference for analysis on symmetric spaces and invari- 
ant differential operators is Ref. |22]. The appendices in 
Ref.|23 contain an accessible introduction to root systems 
and symmetric spaces. 



1. Class BD 

In this case the space of transfer matrices, of type DI 
in Cartan's classification, is M — G/K, where 

G = SO Q (N,N), K = SO(JV) x SO(JV). (Al) 

The Lie algebra of the group G is Lie(G) = so(N, N), 
and it admits the Cartan decomposition: 

so(JV, N) = so(N) + ao(N) + V, (A2) 

where the subspace V consists of matrices of the form 



9 

r o 



(A3) 



with a real N x N matrix 9. 

Within the space V there is a maximal Abelian sub- 
space A, and one convenient choice is the set of matrices 
of the form 



ax 



X 
X 



X = diag(xi, . . . , x N ), x l G 



(A4) 



Let us choose a canonical orthonormal basis {ei . . . e^r} 
in the space A* dual to A (this space is isomorphic to 
M. N as a vector space) so that 

ei(a x )=Xi. (A5) 

Then the set of roots a of the pair (Lie(G), A) is 

Z = {±e i ±e j , l^i<j^N}. (A6) 

These vectors form a root system of type D^y. All the 
roots have the same length, and their multiplicities m a 
are all equal to 



1. 



(A7) 



It is convenient for our purposes to choose as positive 
the roots 



S+ = {e j ±e i! l^i<j^N}. 



(A8) 



This somewhat non-standard choice leads to the follow- 
ing dominant Weyl chamber (note that x\ may be posi- 
tive or negative): 



xi < x 2 < 



< x N . 



(A9) 



This makes the middle element in the Cartan decompo- 
sition Q unique. 

The half-sum of the positive roots (|A8|) is 



1 N 



p(x) = p(a x ) = - l)a 



(A10) 



The density function (usually denoted 8{x) in the math- 
ematical literature^) 



J{x) = J^J (sinha(ox))' 



= J^J sinh(xj + Xi) sinh(iEj — x») (All) 

i<j 

is exactly the Jacobian which appears in the Fokker- 
Planck equation (|24|l . 

The Weyl group W of the root system of type D^r is the 
semidirect product of the symmetric group Sn (permut- 
ing the and Z 2 _1 (acting by an even number of sign 
changes ei — * — e/). The ring of polynomial invariants of 
W is generated by 



S2, Si, . . . , S2N-2,PN, 



(A12) 



where 



:2k 



PN 



(x) = Y[x l . (A13) 



We can also form Vl^-invariant differential operators on 
A replacing Xi by di = d/dxi in these expressions. The 
operators obtained in this way, 



s 2 (<9), . . . , S 2 N-2(d),PN(d) 



(A14) 



are local analogs (valid near Xj = for all i) of the gen- 
erators 



A S2 , . . . , A S2JV _ 2 , A pjv 



(A15) 



of the algebra of A'-radial parts of the globally G- 
invariant differential operators on G/K. This may be 
understood directly, by comparing these expressions with 
Eqs. (|35l I36fl . transformed to radial variables, and drop- 
ping the angular parts. 

In general, only A S2 = A2, the radial part of the 
Laplace-Beltrami operator is known explicitly: 



J2 J 



dxi dxj 



(A16) 



Conjugation of this expression by the density function 
leads to the operator D 2 appearing in the FP equation, 
and the other operators D S2k = D 2 k and D pN are ob- 
tained similarly; see Eqs. (|27|) or l|39|l . 



18 



For the purposes of this paper we do not need explicit 
expressions for the Ds, but only their asymptotic limits in 
the regime of a long wire, when the radial coordinates are 
widely separated, \xi\ <C x% -C ■ ■ . -C xjv- In this regime 
many simplifications occur. First, the density function 
can be approximated as 

«e 2 "W =J| e 2(i-i)x 4 . (A17) 

i 

Secondly, it follows from Theorem II. 5. 23 in Ref. I22I that 
the differential operators of Eq. (|A15|I tend to 

A Sk A Sk = e -'W Si (fl)e'W - e-^)( Sfe (a)e^) 

JV 
i=l 

A PN A PJV = e-^W(S)e^) - e^^^e**)) 

JV JV 

= ]J e-^flfe^ - n 0<p(a:). (A18) 



i=i 



i=l 



In the subtracted terms the derivatives act only on 
e p ( x \ producing x-independent constants so that no non- 
derivative terms appear in the operators. 

Then it follows that the asymptotic forms we need are 

D Sk -» D Sk = e»M Sk (d)e-rM _ e -p(*>( afc (0) e '(*)) 

JV 

= S ((e p(K) ^e- p(a:) )' c - (i - 1) A 

2=1 

- D PN = e P^ PN (d)e-^ - e~P^( PN {d)e^) 

N N 

= Y\e p ^d l e- p{x) -JJ(i-l). (A19) 



For class BD we have (see Eq. IjAlOp ) 

e P(x) d . e - P (x) = g._ i + lj 



(A20) 



so, for example, the FP equation in the asymptotic limit 
contains 

N 

5 2 = ^r((^-*+i) 2 -(z-i) 2 ) 

i=\ 
N 

= Y,(d?-2(i-m). (A21) 



Finally, the operator we are interested in most, D PN , in 
the asymptotic limit is replaced by 



JV 



D PN = l[(d t - i + 1) = di(d 2 -l)...(d N -N+l). 

i=l 

(A22) 



2. Class Dili 

In this case the space of transfer matrices, of type D 
in Cartan's classification, is M — G/K, where 



G = SO(2jV,C), K = SO(2N). 



(A23) 



The Lie algebra Lie(G) = so(2jV, C) of the group G is 
usually taken to consist of all antisymmetric complex 
2N x 2N matrices. The Cartan decomposition 



so(2N, C) = so(2A0 + V 



(A24) 



coincides with the decomposition of antisymmetric com- 
plex matrices into real (for so(2iV)) and imaginary (for 
V) parts. 

By the unitary transformation with 



S 



1 (In Hn 



y/2 \In -Un 



(A25) 



(In is the N x N unit matrix) the antisymmetric complex 
matrices may be replaced by matrices of the form 



C B a t), B 1 = —B, C* = -C, (A26) 



where A, B, C are complex N x N matrices. In this basis 
the Cartan decomposition can be described as follows. 
Elements of so(2jV) are matrices l|A26(l with antisym- 
metric real part and symmetric imaginary part, and vice 
versa for the elements of V . 

A maximal abelian subspace A of V is chosen as the 
set of matrices 



X 

a x = I r, v ], X = diag(xi, . . . ,x N ), x, t £ 



-X 



(A27) 



It follows then that the root system in this case is the 
same as for class BD, Eq. (|A6|) . except that the multi- 
plicities are now 



(A28) 



The positive roots and the dominant Weyl chamber are 
the same as before, so most of the expressions from the 
previous subsection apply here as well. Some changes are 



JV 



p{x)=^22(i-l)x if (A29) 

i=l 

J(x) = J^J smh 2 (xj + Xi) sinh 2 (a;j — xA. (A30) 



The main difference, however, comes from the fact (which 
is true whenever the group G is complex) that when all 
the roots have multiplicity 2, the function J 1 / 2 (a;) is a lin- 
ear combination of exponentials of positive roots. Then 
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it is possible (see Theorem II. 5. 37 in Ref. 22) to find ex- 
plicit expressions for all the radial parts of Eq. (1 A 15(1 : 

A Sfc = J- 1 /2 Sfc(a)J l/2 _ J-l/2( Sfe(9)J V2) 5 

J- l l 2 PN {d)J 1 ' 2 - J-^ipMid)^ 2 ). (A31) 



As before, the subtractions remove constant non- 
derivative terms from the differential operators. 

Correspondingly, the generalized Fokker-Planck equa- 
tion contains 

D Sk = J^s^d)^ 2 - J-y\ Sh {d)J^), 

D PN = J^ 2 PN (d)J-^ 2 - J- l l 2 (p N {d)J 1 ' 2 ). (A32) 

As in the previous case of class BD, we do not really need 
these explicit expressions, and, instead, use the appropri- 
ate asymptotic limits: 

N 



D2 = j2(%-^-m), 



i=l 
N 



D PN = Y[(d t -2i + 2) 



The Fokker-Planck equation (|24(l can be solved in 
general (for any symmetry class) with the help of the 
so called spherical transform on the space G/K, see 
Ref.lU The solution is expressed in terms of the so called 
spherical functions <pk(x) labelled by N real parameters 
k = {ki, . . . , fcjy}- These functions are the eigenfunc- 
tions of A2 with eigenvalue —k 2 = — kf. The solution 
for P(x; L) (up to a L-dependent normalization constant 
N(L)) in terms of the spherical functions is an analog of 
the Fourier transform on the double coset K\G/K: 



P(x;L) = N(L)J(x) 



\c(k)\ 2 



CXP (-2^ fe(a:) - 

(B5) 



The Harish-Chandra function c(k) appearing in the inte- 
gration measure is known in all cases, and for the BdG 
classes is given by 



N 



c{k) = \. 



N 

nn 

i<3 ± 



r( 



(B6) 



d 1 {d 2 -2)...{d N -2N + 2). (A33) 



APPENDIX B: SOLUTION OF THE 
FOKKER-PLANCK EQUATION 

Equation (|24|l is supplemented by the initial condi- 
tion which corresponds to M = 1 for L = 0. There 
are also some boundary conditions imposed on the func- 
tion P(x;L). They follow from the fact that the Fokker- 
Planck equation l|24|l has the form of the continuity equa- 
tion for the probability current Si (we use the short-hand 
notation 8l = d/dL, di = d/dxi): 



d L P 



1 N 



2~{l 



Si = diP-PdAnJ. 



(Bl) 



Then the conservation of the total probability 
J c dxP(x;L), where the integration is over the do- 
minant Weyl chamber C, is ensured if the normal 
components of the current vanish on the boundary of C: 



Sn\ xe9c - °- 



(B2) 



The boundary dC is a union of hyperplanes 2^+1 = Xi 
for i = 1 through N — 1, and also x 2 = —x±. Thus we 
have, for example, 



{d t P - Pdi In J) - (d i+1 P - Pd i+1 In J) 



= 0. 

(B3) 



In addition, as we will see momentarily, the function 
P(x; L) itself also vanishes on dC, see discussion after 
Eq. p6ll: 



P(x;L) 



xedc 



= 0. 



(B4) 



From the solution (|B5|I it follows that the distribution 
density P(x; L) vanishes on the boundaries of the Weyl 
chamber C, because the Jacobian J(x) does, and because 
the spherical functions ifk(x) are bounded (see Theorem 
II.8.1 in Ref. 12). 

The spherical functions <fik(x) are known explicitly 
only for classes A, AIII, CI, and Dili in the AZ clas- 
sification, and so in these cases there are explicit expres- 
sions for P(x; L). For all the symmetry classes, however, 
the expression l|B5|) can be significantly simplified in the 
asymptotic limit of long wires L ^> j£. In this case the 
"diffusing" radial coordinates become very large, ^ > 1, 
and for such values of Xi the spherical functions tend to 
(see Theorem II.5.5 in Ref. HI) 



ip k {x) ~ J 1/2 {x) ^2 c ( sk ) e ' 



(sk,x) 



(B7) 



where the summation is over the elements s of the Weyl 
group W, and (, ) denotes the scalar product. It is im- 
portant that for classes BD and Dili this is still valid 
even if \xi \ is not large. In the limit L ^> ~f£, only ki <§; 1 
significantly contribute to the integral in Eq. (|B5|) . For 
small values of k-, the Harish-Chandra function becomes 



N N 
,S(mi) ■ 



c-\k)ocl[k^l[(k]-k 



I I ■ 



(B8) 



i<3 



where 



r/ n f 0. for mi = 0, 
6 ^ = { I! for m \ > 0. 



Then the integral in Eq. (|B5(1 can be done, and with the 
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use of Eq. i'25\) one can obtain 



P(x; L) «JV(L) f[ exp (- ^x\) xf mi ] (sinh 2x t ) 

i=l 

N 

x - x 2 )(sinh 2 Xj - sinh 2 a;,)" 1 "/ 2 



m ; /2 



(B9) 

Further simplification is possible if we use the fact that 
xi (i > 1) are not only large in this regime, but they are 
also widely separated: 

|xi| < x 2 < . . . < x N . (BIO) 

In this portion of the dominant Weyl chamber far from 
the origin and the chamber's boundaries the asymptotic 
solution 1|B9|) simplifies further to the form given in the 
main text, Eq. 



APPENDIX C: DETAILS OF THE SUSY 
SOLUTION 

The supersymmetry technique is most useful for the 
calculation of moments of Green's functions. The latter 
can be viewed as sums over paths, each path contribut- 
ing a product of factors that are the amplitudes (ele- 
ments of the Hamiltonian or transfer matrix) for each 
step of the path. Hence it is natural to represent these 
sums in terms of Green's functions in a field theory for 
particles (fermions and bosons) that propagate accord- 
ing to these amplitudes. With equal numbers of fermion 
and boson components, the partition function is unity, 
and the averages of products of Green's functions can be 
calculated 28-34 35-36 . The resulting field theories possess 
symmetries (supersymmetry) that rotate fermion and bo- 
son fields into each other. This symmetry guarantees that 
the partition function is indeed unity. 

In the present case, for class BD, the two directions of 
propagation along the wire will be viewed as vertical and 
labeled "up" and "down" respectively. The logarithm of 
the transfer matrix will effectively become a Hamiltonian 
that propagates the particles forwards (up). At each time 
(i.e. position along the wire, as in the main text), the 
two sets of variables, up and down, will be grouped as 
two "sites" . We will view this system then as a two-site 
quantum Hamiltonian (in imaginary time). For technical 
reasons, it is necessary at some stage to treat the cases 
N even or odd separately. We consider only iV even, as 
it describes superconducting wires, and write N — 2n for 
convenience. 

On the up site, we may choose the bosons and fermions 
bai fa (Latin indices run from 1 to N, and we assume 
summation over repeated indices) to obey canonical com- 
mutation relations. On the down site the bosons b a are 
canonical, but the fermions f a satisfy the following com- 
mutation relation: 



As the dagger f always stands for the adjoint, this re- 
lation results in negative square norms for states with 
odd numbers of / (if we choose the vacuum to have pos- 
itive norm). This is necessary for the cancellations that 
result in the partition function being unity (before aver- 
aging over disorder). (Other choices for where the minus 
signs go are possible, but do not essentially change any 
results.) We remark that we could use J\f copies of the 
fermions and bosons (Af = 1, 2, . . . ), and this is neces- 
sary when studying higher moments of Green's function 
than we will need here. 

For the transfer matrices, it is convenient for us here to 
use the Cartan decomposition for the transfer matrices 
at the nodes in the form 



T = kp. 



(C2) 



The first factor here k e K = SO (AT) x SO (AT) describes 
the propagation and mixing, or forward scattering, of 
particle fluxes on the up and down sites (the two factors 
of SO(N)xSO(N) correspond to the up and down sites 
respectively). As usual, if we assume the uniform distri- 
bution for k with respect to the Haar measure on K, then 
averaging over this distribution projects the fermionic 
and bosonic Fock spaces to the subspaces of the singlets 
of SO (AT) on the up and separately on the down sites. 



1. The transfer matrix and the Hamiltonian 

The second factor p 6 expP in the Cartan decompo- 
sition l|C2|l describes the evolution of the states across a 
node of the network. The element (|14|l of V parameter- 
ized by 9 gives rise to the transfer matrix 



6 
ex P I aT o 



F sinh Z^ T coshv^ 



(C3) 



where the matrices 9 T 9 and 99 T have non-negative eigen- 
values, and the square roots are well defined. 

The second quantized version of the transfer matrix p 
is the evolution operator 

V = exp (i6 ab J ab ), (C4) 
J ab = blb\ + hb a + fl ft + f b f a . (C5) 

As before, we assume that the transfer matrices are 
random with some generic A'-invariant distribution P(0). 
As in our analysis of the FP equation, it will be sufficient 
to consider only the cumulants and l|20l) . Then, 
averaging V over the distribution of the disorder P{9) 
we write it as 



{fa> ft } — 



-Sab 



(CI) 



[V] = 



-aH 



(C6) 
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where the effective Hamiltonian H contains products of 
J a b with cumulants of 9, and a is the time step, as before. 
The two non-trivial cumulants Ijlfil 121 H give two terms in 
the Hamiltonian (we choose to rescale the Hamiltonian 
by a factor): 



H r — - — H — H + 2XHi, 
2c 2 



(C7) 



Hq — -^5 ai a 2 5b 1 b 2 Ja 1 biJa 2 b2 — ~^JabJabi (C8) 



Hi 



N i 



N 



2 (TV!) 2 
where the parameter 



b N Ja-^bx ■ ■ • Ja N b N j (C9) 



(CIO) 



is chosen to simplify equations in the following. Note that 
this dimensionless parameter is related to physical length 
scales introduced in Section ITTT1 as |A| = £x/(2£typ)- 

As usual, we should add to the Hamiltonian H r a reg- 
ulating term 



"(bX + flf a + b\b a - /t/J, 



(Cll) 



which suppresses large boson numbers, to represents the 
imaginary part of the energy e in the original random 
Hamiltonian, and is necessary when studying e.g. the 
DOS. The full effective (rescaled) Hamiltonian is then 



H r = Hq + 2XHi + H u 



(C12) 



2. Supersymmetry algebra osp(2|2), and its 
representations R and R 

The subspaces of states that are invariant under 
SO(-/V)xSO(iV) are infinite-dimensional representations 
of the superalgebra osp(2|2), which we will see is the 
symmetry algebra of the problem. We denote by R the 
representation space on the up site, and by R the one on 
the down site. 

The osp(2|2) symmetry is generated by the bilinears 
in the fermions and bosons that are invariants of SO (TV). 
The generators of osp(2|2) on the up site (in the represen- 
tation R) are given by (using notation for the generators 
from Ref . y2 



B = 


\ (flfa - n ) , 


Q 3 - 


\ (b\b a 


+ n) , 


Q+ = 


l(bl)\ 


Q- = 


-Ql = 


-\bl 


v+ = 


V2 Ja Q ' 


W- = 


-vt = 


-T2 bJ( 




V2 Ja °' 


W+ = 


-vl = 


V2 aJa 



Similarly, on the down site (in the representation R) the 
generators are 



B = 



V+ = 



V- = 



2 V J a J a 1 

\bl 

V2 a/a ' 



Q- 

w+ = - 



-\(b\b a +n)> 

.yt - Lfct ft 

_ \[2 

V2 Ja a 



(C14) 



The generalization to osp(2A/"|2jV) is straightforward, 
but not needed here. 

Using these generators we now describe the representa- 
tions R and R. The representation R splits into two irre- 
ducible representations of the algebra osp(2|2). The first 
one of these, which we denote as Rq is a lowest weight 
irreducible representation with the lowest weight state 
being the vacuum |0) for all b a ,f a . Other states are ob- 
tained by repeated actions of Q + and V±. Normalized 
states obtained in this way are 



\2m,0) = a m Q™\0), 
|2to-1,1) = b m V^\2m,0), 



where 



(n-1)! 



to! (to 



D! 



1/2 



= m- 1 ' 2 . 



(C15) 



(C16) 



The other subrepresentation of R denoted Rm is also a 
lowest weight irreducible representation of osp(2|2), but 
now the lowest weight state is 



|0,7V) = / 1 t / 2 t .../t r |0). 

Other normalized states in Rn are 

\2m,N) = a m Q™\0,N), 
|2m+l,JV-l) = c m W + \2m,N), 

with a m as in (|C16|) and 

c m = {m+ n)~ 1/2 . 



(C17) 



(C18) 



(C19) 



The representation R is the dual of R. Again, it is split 
into Rq and Rn under osp(2|2). Both these subrepresen- 
tations are highest weight irreducible representations. In 
-Ro the highest weight state is the vacuum |0) for b a ,f a , 
and the other normalized states are 



2to 



2m, 0) =a m Q™\0), 
M) = b m W + \2rr^0), 



(C20) 



with the same a m ,b m as in (IC16I). Note that the states 



with odd number of fermions f a have negative squared 
norms. For example, 



(2to-1,1|2to-1,1) = -1. 



(C21) 
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q I2m,0> q + 



I0> 
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I2m-l,l> 



-N/4 



I2m+1,N-1> 




V. 



Q- l2m,N> Q+ 



IO,N> 



IO,N> 



e. l2m ' N> q + 




I2m+1,N-1> 



a 



I2m-l,l> 




io> 



I2m,0> y+ 



FIG. 2: Weight diagram for states on the representations R and _R. 
In i?7v the highest weight state is 



site. In the representation R we have 
1 



and the rest is 



|o,jv> = /7jg.../iio) 



2m, N) =a m Q m \0,N) 7 



(C22) 



|2m+l,iV-l) = c m V-\2m,N). 



(C23) 



The states |2m+l,7V— l) are negative norm states. 

The osp(2|2) weights of the states in R and R are 
shown in Fig. The negative norm states are shown 
with filled circles. We also show the action of the osp(2|2) 
generators on the states. We note here that the bosonic 
part of the supergroup OSp(2|2) under consideration is 
0(2)xSp(2, R), and 0(2) has two connected components. 
In the same way that 0(2) has two-dimensional irre- 
ducible representations [whereas SO (2), and its Lie al- 
gebra so(2), have one-dimensional irreducibles] , the two 
subspaces Rq and Rn together form an irreducible repre- 
sentation of OSp(2|2) 29 . The same is true for the states 
in R. 

The term J^ b in H is an osp(2|2) invariant, and is pro- 
portional to the osp(2|2)-invariant product of the gener- 
ators acting in the representations R and R. It is conve- 
nient to choose the constant c 2 = 1/2, and then 



H = 2BB - 2Q 3 Q 3 - Q + Q_ - Q_Q+ 

-V+W- + W-V+ +V-W+ -W+V-. (C24) 

Then it is clear that Ho acts irreducibly within each of 
the four subspaces Rq <g> -Ro, ■ ■ ■ , Rn <8 Rn- 

The term Hi has a more complicated structure. It 
is an osp(2|2)-invariant product of N operators on each 



1 

M' 
l 

M' 
l 

1 



ft ft 


C- 




f ft ft 


D- 


= Dl 


b ft ft 


R- 


= R\ 


b f f 1 f f 


Y- 


= 4 


b ft ft 
aw "ai J a 2 "'J a N ' 


Y+ 


= xi 



In the representation R we have 
C = — — f ft ft 

w — c ai ...ajy J a\ " 'J ajy ) 

D = ——F f ft ft 

= ^j £ ai...a N ^ 1 ^ 2 /] 3 "-/i s i -R+ — -R- 

Y = ——F P ft ft 

± — ^yj G «l ...ajv "ai J 02 " 'J a» ' 

y, = — f ft ft ft 



C+ = ci 



X + = Y 1 



X-=Y1. (C26) 



Then Hi projected onto R ® i? is proportional to the 
osp(2|2) invariant product of these operators, given by 

Hx = — [C+C- + C-C+ - N(D + D- + D-D+) 

- N(N - l)(R + R- + R-R+) 

-N{X+Y- - Y-X+ +X-Y+ - Y+X-)]. (C27) 

The terms of Hi that have non- vanishing matrix elements 
within R<E>R flip the states between the subspaces Ro^Ro 
and Rm <g> Rn- 
The term 



H u = 2lu(B + Q 3 -B-Q 3 ) 



(C28) 
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breaks the supersymmetry from osp(2|2) down to a u(l|l) 
subalgebra, which is generated by 



Q 3 ,B,V.,W+. 



3. Eigenvalue problem 



(C29) 



For non-zero to in the symmetry of the Hamiltonian 
H r in Eq. is u(l|l) of Eq. (|U29)l . We can use this to 

find eigenstates of H r in the singlet sector of this u(l|l). 
Basis states in this sectors |m)j are defined in terms of 
the states in R <8> R as follows: 



|m)i 



m > 1, 



|0)|0), 

-^=(|2m,0>|2m70) 
+ |2m-l, l>|2m-l,l>V 

-*0 

+ |2m-l,iV-l)|2m-l,A/-l}Y m > 1. (C30) 



We look for eigenstates of if of the form 



|m; 2 = — I \2m-2, N)\2m-2, N) 



.328ft , we obtain the action of the parts of the Hamilto- 
nian H r . 



flo|0}i = -nV2|l)i, 
Z?o|to)i = — (m + l)(m + n)|m + l)i 
+ m(2m + 2n — 1)|to)i 

— (m — l)(m + ji — l)|m — to ^ 1, 
Ho\m)2 = —m(m + n)\m + 1)2 

+ (2m — l)(m + n — 1)|to)2 

— (m — 1)(to + n — 2)|m — 1)2, to 1. 

(C32) 

fli|0)i=m/2|l) 2 , 

1 77T.) 1 = (to + n)|m + 1)2 

— (m +n — 1)|to) 2 , to 1, 
1 777.) 2 = to|to)i — (m — 1)|to — 1)i, to ^ 1. 

(C33) 

H u \m)i = 2tow|to)i, 

H u \rn)2 = 2(m + n — l)ui\rri)2, Vto. 

(C34) 

Note, that in the subspace spanned by the states |m)j 
we have 



\<l>) = ih\oh 



E^Ni + ^Wa). ( C31 ) 



Hl\m)i = -H \m)i. 



(C35) 



where the coefficients V'Oj V'm j and V% satisfy some re- 
cursion relations, and can be found from them. The re- 
lations are obtained by plugging the ansatz l|C31|) into 
the eigenvalue problem H r \ip) = E r \ip) and equating the 
coefficients of basis states on both sides of the equation. 

First, using the definitions (IHHl IUl4l IH251 of 
the operators appearing in the Hamiltonians (|C24I IC27I 



(It is to have the simplest possible coefficient in this equa- 
tion that the constant A was introduced as in Eq. I|C10|) .) 
This relation will significantly simplify the analysis of the 
recursion relations for ip m for small to (see section Rj 51 . 

Combining the above equations, we get the following 
recursion relations: 



E r lpO = 0, 

{2w-E- R)ip[ 1] = (n + l)^ - (2n + + nv^o + 2A(^ 2) - ^i 2) ), 

(2ncj - Er)^ = n(4 2) - V>i 2) + 2A(V4 1} - V^o)), 
-E r (A^ + = (to + 1) ((to + n + 1)A 2 V^ } + (1 - 2w)A^&> - 2u^&> + 2A(A 2 ^ } + A^ (2) )) , "i ^ 2, 

-£ r (A^ + ^) = (to + n) ((to + 1)AVL 2) + (1 - 2co)A^ ) - 2^ 3 + 2AA^ } ) , to ^ 2. 

(C36) 



The last two equations are written in terms of the finite 
differences: 

Alp m = 1p m+ l ~ 1p m , 
A 2 1p m = lpm+2 - 24>m+l + 1pm, (C37) 



Our strategy for the analysis of these equations is the 
following. For large m 3> n we can treat to as a continu- 
ous variable. The solutions of the resulting equation are 
required to decay at infinity. The solutions of the contin- 
uous equation show that the effect of small u> becomes sig- 
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nificant only for m > uj^ 1 . Then for m <C uj^ 1 we can ne- 
glect oj in the equation HC36|) . The corresponding Hamil- 
tonian (within the space R ® R) is simply 2XHi — B.\ . 
Then, instead of solving Eq. (|C36(I . we can first find the 
eigenstates of H\ and then combine them to obtain the 
energy eigenstates. 

After finding approximate solutions in two regions 
m n and m <C we can match them in the 

overlapping region n <C m <C lu^ 1 (we choose u> to be 
asymptotically small so that this region is large). For 
the ground state this matching procedure gives a com- 
plete solution. We can use this solution to calculate the 
expectation value of the "order parameter" in the ground 
state, which gives the density of states in the ID wire. 

For excited states the matching procedure gives the 
spectrum of energies E of the coupled superspins. Then 
the localization length £ mca n associated with average cor- 
relators in the state at the Fermi energy in the original 
ID wire problem is given by the inverse of the gap in the 
energy spectrum above the ground state in the superspin 
problem. Let us implement this strategy step by step. 



4. Solution of Eqs. (IC36I) for large m 

For large m n, m can be treated as a continuous 
variable, and we can replace finite differences by deriva- 
tives with respect to m: 

(C38) 

In this limit we neglect uj and constants of order 1 (in- 
cluding n) compared to m. Assuming a smooth variation 
of ^i{m) we also neglect ifj"(m) compared to ip'^m), and 
ipKm) compared to ipi(m). Then the Eqs. I|C36|) become 

(2a;m — E r )ipi(m) — m(m%pi(m) + ^[(m) + 2\ip' 2 ( m )) , 

(2ojiti — E r )ip2{?n) = m(mip2 (m) + ip 2 ( m ) + 2A?/> 1 (m)). 

(C39) 



where A± are constants, K the modified Bessel function, 
and 



To solve the equtions (|C39ll . we introduce 

i>±{m) = V>i(m)±^ 2 (m). (C40) 

Then, Eqs. (|C39(1 can be written as 

(2wm - E r )ip±(m) = m 2 ip±(m) + (1 ± 2X)mip'±(m). 

(C41) 

After a change of variable x 2 = 8umi, and the substi- 
tution ip±(x) = x T2X f±(x), we get the modified Bessel 
equation 

x 2 fi{x) + xf' ± {x) - (x 2 + 4(A 2 - E r )) f±(x) = 0. 

(C42) 

The decaying solutions for ij)± can now be written as 

2J K 2 ^x), ( C43 ) 



VA 2 - E r 



(C44) 



For the purpose of matching these solutions to the ones 
for m <C we need the small x behavior of the Bessel 
function Kn(x). To obtain it we use the definition 



KM 



-(I^(x)-W) 



2 sin fin 

and the small x asymptotics for I^(x): 
I (x) ~ {X,2Y 

lfi[x) r(i + M )- 

Then we get for small x 



(C45) 



(C46) 



(C47) 



After the substitution x = 2(2u;m) 1 / 2 the limiting form 
of tp± for small m becomes 

V>±(m) ~ A± [T(2/j l )(2ujm)~^ x + r(-2^)(2mw)^ TA ] . 

(C48) 

For the ground state E r — 0, we have fi = |A|, and the 
asymptotic solution 1C48I) takes the form 



ip±(m) ~ A± [r(T2A) + r(±2A)(2mw) T2A ] . (C49) 

Let us specify the previous solutions to the case of 
A = 0. Eq. I|C48|) for an excited state wave function 
simply becomes 

V>±(m) ~ A ± [r(2/i)(2wm) _ " + r(-2/i)(2 

(C50) 

where now /1 = iy/E r . For the ground state we have from 
Eq. iC43l 

ip±(m) = 2A ± K Q (x) A±(ln(2wm) + 2 7 ), (C51) 

where 7 is the Euler's constant. 



5. Solution of Eqs. (IC36I) for small m 

In the limit m <C i-j -1 , can be neglected and thus 
we have 

[2Aff x - Hi] |</>) = EM), (C52) 
where we used Eq. (|C35|) . Given the eigenstates of H\\ 

H 1 \iP,j)=M,j), (C53) 
the eigenvalues j are related to E r by 

f - 2Xj + E r = 0, (C54) 
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which has two solutions 

ji2 =\±(1, 



(C55) 



with the same /i as before, see Eq. l|C44ll . The energy 
eigenstates are then given by 

\i>)=A 1 \Tp,j 1 )+A 2 \i>,j 2 ). (C56) 

Equation (|C53J) leads to the recursion relations 

jipo = 0, 
j4 2) =n(^-Vi 1) ), 

j^W = {m + n-l){^_ 1 -i>$). m^2. (C57) 



First, we consider excited states with E r > 0. Then 
j 0, and it follows from the first of Eqs. (|C57(1 that 



4>o = 0. 



(C58) 



The rest of Eqs. (|C57(1 can be solved using the generating 
functions 



/<(*) = E ^« m - x . 



The coefficients are obtained from the generating 
functions as 



(C59) 



« 



J 27TI 



(C60) 



where the integral is taken counterclockwise around zero. 

From Eqs. I|C57|I . we obtain the equations that the 
generating functions satisfy: 

z(l -z)^ + [n-(l + n)z] fx + jf 2 = 

(l-z)f^-/ 2 +j/i=0. (C61) 
dz 



Substitution of fi(z) from the second Eq. I|C61|I into 
the first one, and the change f 2 (z) = (1 — 2)- ,_1 it(z), lead 
to a hypergeometric equation 



d^u du 
z(l - z)— + (n - (n + 2j + l)z) — - j(j + n)u = 0, 



dz 



(C62) 



with the solution which is regular near z = 



u(z) = F(j,j + n;n;z). (C63) 
Thus, in terms of the hypergeometric function 

=A^z- m (l-zy- 1 F(j,j + n;n-,z), (C64) 

where the constant A is determined later by matching 
this solution to the one for large m. 

To evaluate the last integral, we use one of the Rum- 
mer's relations for the hypergeometric function: 



r(n)V(—2i) 

(1 - zy- l F{j, j + n; n- z) = 1 ' \ r( J> .M - zy^Ffj, 1 - n + j; 1 + 2j; 1 - z' 1 ) + (j <-> -j). (C65) 

r(n- j)r(-j) 

Then 

(C66) 



Vr(n-j)r(-j) y 



= ^ - 1 - n + j; 1 + 2j; 1 - z" 1 ). (C67) 

In the last integral we deform the contour to run along the branch cut [1, oo). On the upper and lower sides of the cut 
1 — z = e~ M (z — 1) and 1 — z = e t7r (z — 1), correspondingly. Then the sum of the contributions from the two sides is 

I(j) = E±ll / dzz-t m+: >\z-iy- 1 F{j.l-n + j;l + 2j-A-z- 1 ). (C68) 

T Jl 

The change of variables x = 1 — z~ x leads to a tabulated integral: 

1 r(i)r(m) 

dx x'-\l - xr^FU, 1 - n + j; 1 + 2j; x) = { ' ) ' 3 F 2 (j, 1 - n + j, j; 1 + 2j, m + j; 1), (C69) 

r(m + j) 

I 



see Eq. 7.512.5 in Gradshteyn and Ryzhik2£. This pro- vides us with the exact solution of the recursion relations 
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j571 . However, we really need only the asymptotic form 
of this solution for large m: 



m 

Denoting 



1 



r(i - j 
Mj) 



-m~ j (l + Oil/m)). 
T{n)T{2j) 



T(n + j)T(j)T(l+jY 

(2) 



(C70) 



(C71) 



we have the asymptotic behavior of ipm for large m: 

~ A(k[j)m? + A(-i)m^). (C72) 

Next we substitute this asymptotic form into the recur- 
sion relation for ipm in Eq. I)C57|1 . and obtain with the 
same accuracy 

r m 



-A(A(j)m j - A(-j)m j ). 



(C73) 



For a given j we get the asymptotics for the functions 



tp±(m) ~ ±2AA(Tj)m Tj . 



(C74) 



Combining these for both values of j, we get, finally, for 
excited states 

^±(m) ~ ±2(A 1 A(TJi)m Tn + A 2 A(T32)m T32 )- 

(C75) 

Next, we turn to the ground state. In this case E r = 0, 
and we have 



h = 2A, 32 = 0. 

For j = the Eqs. I|C57|) give 
tpo = const, 

vi 2) =4 2) . 

Then 

V>±(m) = V2^ ±4 2) - 



(C76) 



(C77) 
(C78) 



Combining this with Eq. (|C74J) we find that the in the 
case A ^ the Eq. I|C75|I specializes to 

V»±(m) ~ ±2AiA(t2A)to =f2A + V2ip ± V^- (C79) 

In the special case of A = the roots ji and 32 become 
degenerate, and we have to be more careful. In this case 
it is actually easier to come back to the recursion rela- 
tions l|C36|l . The first of these leaves ipo an unspecified 
constant, while the rest can now be written as 



(n + 




= n 




4 2) - 4 2) = 0, 


A 




= 0, 




A mA^ 1 


= 0, 



m Ss 2, 
m > 2. 



(C80) 



Integrating the last pair of equations i|C80|l . and finding 
integration constants from the first pair, we obtain 



(m + n) 



in 



\b {1) 

W,( 2 ) _ ^(2) 
V^m+l 1 J m 



= n(4 1] - V24>o), 

= 0. (C81) 



Thus, 



n + m — 1 

= \/2> + n^i 15 - v^V-o) + m) - *(n)] , 

(C82) 

(2) 

where ^ is the Digamma function. For ipm on the other 
hand, we have 



",(2) 



(C83) 



vi 2) = < 

For the combinations ip± this gives 

ip±(m) = V2V-o±4 2) 

+ - \/2> ) [*(n + m) - , (C84) 

with the limiting form, for large m, 

+ - V2ip ) [him - *(n)] . (C85) 

6. Matching solutions for large and small m 

Let us start with the ground state E r — and a generic 
A 7^ 0. In this case we need to match the Eqs. (|C491 and 
(|C79|I . This gives the conditions 



A±r( T 2A) = \/2^o±V>i 2) , 
A±r(±2A)(2w) T2A = ±2,4iA(=f2A). (C86) 

Solving for A + and A_ , we get 

2V2^o 



A, 



where 



$(A) 



r(2A)(l-$(A)(2w) 4A )' 

272^0 $(A)(2u) 4A 
~r(-2A) l-$(A)(2w) 4A ' 

T 2 (-2A)A(-2A) 



r 2 (2A)A(2A) 



(C87) 



(C88) 



For finite A the leading behavior for asymptotically small 
lu is then given by 



A + cx oj 



4 A 



A_ = const, if A > 0, 



A+ = const, A- cx w 4|A| , if A < 0. (C89) 
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FIG. 3: Solid and dashed lines show the right and left hand 
sides of Eq. l|C96jl for n = 5, A = 0.2, and cu = 0.001. 



For A = 0, matching Eq. (|C85(I and Eq. (|C51[) gives 



A A 



-V2> 
ln(Ccj) [ 



Vw, 

C = 2 e 2 T+*(™). 



(C90) 



Essentially the same result is obtained by taking the limit 
of A — > in Eqs. (|C87|I . Indeed, in this limit 



r(2A) 



lim 



1 

2A' 
1 - {2oj) 4X 
2A 



$(A) - 1, 
= -21n(2w), 



(C91) 



and Eq. (jC90|) follows (but with C — 2 now). 

The final form of the ground state (with the choice of 
^0 = 1) for general A ^ is 



tp±(m) = A±(t(=f2A) + T(±2A)(2r 



= 2A±(2Lum) TX K 2 \(2\ / 2ujm), 

For A = we have, instead, 
t/>±(m) = V2( 1 + 

2y/2 



ln(Cw) 



m(Cu;) 
K (2V2ujm) 



m <C — , 
a; 

m 3> n. 
(C92) 



m <C — 

m 3> u. 
(C93) 



Next we discuss excited states with 
E r = X 2 -fi 2 , 



(C94) 



see Eq. (IC44II . The value of fi may be obtained by match- 
ing the large m and small m limits of the -0 l 's in the 
excited states, Eqs. I|C48|) and i|C75|) and by requiring 
a consistent, non-trivial solution. Unlike in the ground 
state, we now have ?/>o = 0, and thus the system of equa- 
tions for coefficients A±, A\ 2 obtained from the matching 



is homogenous: 

A+T{2^){2uj)-^- x = 2A 1 A(-fi - A), 
A+r(-2 M )(2c^- A = 2A 2 A(/x - A), 
A^T(2 f i)(2co)-^ x = -2A 2 A{-^ + A), 
A-T(2fi)(2ujy i+X = -2AiA(/x + A). 



(C95) 



A non-trivial solution of this system exists only when the 
corresponding determinant vanishes, and this gives 

= JJt 0i - X ] M ^ + X) ,,- (C96) 



r 2 (-2/x)A(-/*- A)A(-/z + A)' 

Before analyzing this equation, let us note that, accord- 
ing to the definition in Eq. (|C44(1 . if /i is real, then it 
must lie in the range 

< n < |A|. (C97) 

It can also take purely imaginary values. 

A quick look at the right hand side (RHS) of the 
Eq. (|C96(I shows that the value |A| = 1/4 is special, and 
we will restrict ourselves to the values |A| < 1/4. In 
this case for the allowed real values of /i, the RHS of the 
Eq. I|C96|) monotonically decreases from 1 to 



> 0. 



r(4|A|)r(n-2|A|) 
r(-4|A|)r(n + 2|A|) 
Then for any n there is a value 



1 



r(4|A|)r(n-2|A|) \ 1/4|A 



" 1 " 2 I r(-4|A|)r(n- 



•2 A 



such that for uj < w„ there are no real solution for ji, 
apart from /i = 0. A typical situation is shown in Fig. [3] 
The absence of real solutions for /i means that the en- 
ergy spectrum of the reduced Hamiltonian H r has a gap 
A 2 above the ground state energy (which is by super- 
symmetry). This, in turn, implies that the mean local- 
ization length in the original disordered wire problem is 
finite 

-l 



= 2[^A 2 



(C98) 



except at A = 0, where it diverges. After substitution of 
A from Eq. I|C10|) . and c 2 — a/(j£), the last expression 
reduces exactly to the one given by Eqs. (|42H44I l4*Sjl . 

For purely imaginary /i, however, Eq. 1C96I) has in- 
finitely many solutions, since (2w) 2 ^ = (2w) 2/i e 27rl ', for 
any integer I. Upon expanding both sides of Eq. I|C96|) 
for small \/j,\, one finds that in the limit of u> <C 1, fi takes 
on a discrete set of values, 



1 = 0,1,2,. 



'2|ln(2w)|' 
which gives for the energy eigenvalues 



E r = A 2 



7T 



(C99) 



(C100) 



4 In (2w) 

For oj — > 0, the energy levels above the gap merge into a 
continuum. 
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7. Density of states 



Expressions ilC92IC93|) for the ground state wave func- 
tion allow us to determine the leading behavior of the av- 
erage density of states v(e) in the wire for asymptotically 
low energies e. The density of states (DOS) u(e) for the 
original problem of a ID superconductor is related to the 
"order parameter" of the superspin problem. The "order 
parameter" is the measure of the amount of the symme- 
try breaking induced by oj in the ground state, and is 
expressed as the ground state expectation value of, say, 
the staggered fermionic number 

6 = B- B= -(2n - n f - nj). 

The density of states is found by analytical continuation 
lj — > ie as follows: 



Ke)cxRc(jO)oL_ i£ J, (C101) 
Given the form of the singlets in Eq. 1|C30|I , we have 
6|0) 1 = n|0) 1 , 

6|m)i = ^(n|2m,0)|2m70) 
V2 ^ 

+ (n- l)|2m-l,l)|2m-l,l)V 
0\m)2 = 



- (rc- l)|2m-l, N-l)|2m-l,iV-l)). 

(C102) 

Using the fact that the states 2m— l,l) and 
|2m— 1,N— l) have negative square norms, we obtain 

(O)o = 4 + ^[(tf-(e 2 " 



(C103) 



Consider now the case A = 0. The sum in Eq. I|C103|) 
can be broken up into two parts corresponding to the 
two regimes, separated by some mo such that n <§; mo <C 
l/u). The sum over m ^ mo gives a finite contribution 
as oj — > 0. The rest of the sum in Eq. I|C103I) dominates 
for small lu, and can be approximated by the integral 



A+A_ / dmK%(2V2r 



mo 
2 



2^ 



win (Cw) 



dx x Kq (x) oc 



win (Cuj) 



(C104) 



Thus, for asymptotically small energy the density of 
states behaves as 



v{e) ~ Re- 



1 



(C105) 



ie\v?{Cie) e|ln 3 e|' 
which is precisely the Dyson's singularity of Eq. J2J ■ 

However, in the case of A ^ 0, we have similarly to 
Eq. (|C 104(1 that the contribution from m > mo to the 
"order parameter" is 



(O) ~ A+A- I dm K% x (2y/2um) ~ u 



4|A|-1 



(C106) 



This contribution dominates the finite contribution from 
m < mo only if |A| ^ 1/4 (for larger values of |A| our 
method is not precise enough to determine the DOS). In 
this range we then get the following asymptotic form for 
the DOS: 



MM 



(C107) 



This is precisely the power-law with a non-universal ex- 
ponent depending on the closeness to criticality, which 
was found for a generic ID system in classes BD and 
Dili by Motrunich et. af in a real space RG study. 
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